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AUSTRACI 


The exiraclioti of Aionratics ft om peiroleum naphtha is an inipoitanl coitntieicial 
operation. Because of inct easing cost of energy and ever-increasing need for pure 
aromatic hydrocarbons, the efficient purification and separation of aromatics is 
considered essential for modem refineries. The solvents used for extraction of aromatics 
now a day’s either have good selectivity with simultaneously low solvent capacity or vice 
versa. Therefore it is valuable to search for new solvents, and to estimate the performance 
of existing solvents. 

The selectivity and solvent capacity are the two primary parameters used for 
comparison of solvents. A peifonnance index based on the above two pmameteis has 
been used to categorize the solvents for aromatic extiaclion. fhc selectivity and srjlveni 
capacity are generally estimated using activity coefficient models. These models require 
proper binary interaction parameters for accurate representation of IdJI These 
parameters are not yet available for all aromatic extraction systems. Furthermore the 
parameters available for almost all the systems are found to be specific to that system, 
and can not be used outside that system. For estimation of parameters for thermodynamic 
models such as UNIQUAC and UNIFAC the ternary LLE data have been used. The 
parameter estimation procedure used here has been modified to estimate the parameters 
simultaneously for different systems involving common pairs. The regressed parameters 
using this technique have been shown to be applicable for other ternary and multi- 
component systems. The regressed parameters then have been used to estimate the 
performance of different solvents for aromatic extraction 
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Chapter 1 

Introduction 

Individual aromatic hydrocarbons such as benzene, toluene and xylenes find 
wide use and are considered essential in the modem organochemical industry. Aromatic 
hydrocarbons are key petrochemical raw materials and their importance and versatility 
seem to increase with time. From the light aromatics which are mainly benzene, 
toluene, xylene and generally termed as BTX, benzene is one relatively bright spot in 
the prevailing gloom of petrochemical industry. Since 1960 styrene, phenol and 
cyclohexane have clearly dominated the markets for benzene. Anthraquinone, DDT, 
detergent alkylate, solvent, maleic anliydride and aniline have vied for a position in this 
mai'ket. The other important uses of benzene today are nylon, Benzene-hexachloride, 
2,4-D, bisphenol, dichlorobenzenes, pentachlorophenol and aspirin production. World 
benzene demand is around 31.7 MT and increasing at the rate of 3.6% per annum. In 
India, total supply is about 0.31 MT while its demand is about 0.85 MT. 

Other aromatic compound toluene, once thought of largely as the basis for TNT, 
is now finding expanding use as a solvent and as a chemical intermediate in the 
production of diisocynates, benzyl chloride, methyl styrene, saccharin and terphthalic 
acid etc. Toluene is also used as a blending stock for gasoline. The third aromatic 
compound xylenes also have many uses. Oxidation of various xylenes to phthalic, 
isophthalic and terephthalic acids is typical of new and growing uses of aromatics. The 
other uses of aromatics are dehydrogenation of ethylbenzene to produce styrene 
monomer. 

The major sources for these aromatics from petroleum include straight -run or 
hydrogenated pyrolysis gasoline and catalytic reformats. These sources, apart from 
aromatics, contain large amount of other components also. The compound covers wide 
range of molecular weights and families (Paraffins, Isoparaffins, Olefens, Naphthenes, 
aromatics). Recovery and purification of these aromatic hydrocarbons is therefore an 
important commercial process. It is impossible to obtain a high purity of aromatics by 
distillation, which is the most widely used separation process in petroleum industry, 
due to overlapping boiling ranges of feed components. Many homogeneous binary 
azeotropes exist between aliphatic and aromatic hydrocarbons (Yorulmaz and 
Karpuzcu, 1985). Therefore liquid-liquid extraction will be a better choice 
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provided that components of the mixture show a suitable absolute solubility and 
sufficient difference in solubility in a solvent. And indeed the need to separate mixtures 
of aliphatic and aromatic hydrocarbons provided one of the first large-scale 
applications for solvent extraction. The original Edeleanu process used liquid sulfur 
dioxide for the removal of aromatics from lamp oil. The first major advance in aromatic 
separation occurred with the development of UDEX process, which used mixtures of 
polyethylene glycol and water as the solvent. Processes based on solvents Sulfolane, N- 
Methylpyrrolidone, dimethyl-sulfoxide and N-formylmorpholine followed this process. 

The economics of solvent extraction process largely depends on the selection of 
solvent. The solvents in use now a days for aromatic extraction either have good 
selectivity for aromatics with simultaneously low solvent capacity or vise versa. 
Therefore it is customary to search for new solvents as well as to evaluate the 
performance of existing solvents. To evaluate the performance of solvents the 
important parameters are solvent selectivity, solvent capacity and solvent loss. 
Selectivity signifies the relative separation of two components, while the solvent 
capacity and solvent loss signifies the required flow rate of solvent and the required 
make up solvent respectively. By definition these parameters are calculated using the 
activity coefficients. The liquid phase activity coefficients are calculated using the 
thermodynamic model such as UNIQUAC, UNIFAC or NRTL that can give accurate 
determination of activity coefficients. These models require proper binary interaction 
parameters which can represent liquid - liquid equilibria for highly non - ideal liquid 
mixtures which are normally encountered in aromatic extraction. 

Apart from evaluating the performance of solvents the methodologies known as 
molecular design of solvents, that are used for development of new solvents also 
require the accurate prediction of liquid - liquid equilibria. Therefore the availability of 
binary interaction parameters that can predict LEE accurately is of paramount 
importance for solvent selection. These parameters are yet not entirely available for 
systems encountered in aromatic extraction. For simulation and design of extractor also 
the availability of these parameters is a compulsion. Process simulators like ASPEN 
PLUS may produce drastic errors if used without accurately knowing these parameters. 

These binary interaction parameters are generally estimated using experimental 
liquid - liquid equilibrium data. In case no experimental LLE data for systems of 
concern are available, the infinite dilution activity coefficients can be used for 
parameter estimation, although at the cost of accuracy. The methods used for parameter 
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estimation are least square minimization and maximum likelihood approach. In the 

former technique the weighted errors between experimental and predicted compositions 
is minimized while in the later approach a likelihood function is maximized. The 
drawbacks of least square minimization is that the proper weights are not known at the 
start of regression and improper selection of weights leads to unreliable parameters. 
The maximum likelihood approach, used nowadays, also has a serious drawback. In 
this approach the variances are fixed, generally on the basis of experimental precision. 
But the final values of variances after regression have been found to be inconsistent 
with fixed variances. Therefore this method may also lead to unreliable parameters. 

It is possible to overcome this problem of maximum likelihood approach by 
recomputing the variance at each iteration of the optimization procedure. This 
recomputation of variance automatically weights the objective function thus promising 
most appropriate values of parameters are obtained. 

It has been observed that the parameters obtained using different initial guess of 
parameter lead to different values (Vasquez and Whiting, 2000 [79] ). This is generally 
due to several local minima’s in the objective function used for parameter estimation. 
Furthermore the binary parameters obtained for the same pair are found to be different 
for different ternary systems (Sorensen and Arlt, 1979 [72] ). For example the 
interaction parameters between pair Hexane - Benzene in system Hexane - benzene - 
TMS are different than in the system Hexane - Benzene - TREG as reported by 
(Sorensen and Arlt, 1979 [71] ). Therefore, these parameters are specific to the ternary 
system from which they have been estimated. And hence these parameters can not be 
used for other systems or for the multicomponent LLE prediction. 

'fo avoid these problems, the parameter estimation approach is modified. Here 
the interaction parameters for different systems involving common pairs are estimated 
simultaneously. The objective function is modified to estimate the parameters 
simultaneously for different systems. In this way the problem of getting different sets of 
parameters for same pair has been eliminated. The regression approach is modified so 
that the parameters estimated from ternary LLE data are applicable for multicomponent 
prediction also. 
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Objectives of this Thesis 

1. To establish the methodology for UNIQUAC binary interaction parameter 

estimation from binary and ternary Liquid - Liquid Equilibrium data so that the 
parameters regressed should be applicable for multi-component LLE prediction (for 
aromatic extraction systems). 

2. To establish the methodology for estimation of LTNIFAC group interaction 

parameters from binary or ternary LLE data so that the parameters regressed can be 
used for multi-component prediction (for aromatic extraction systems). 

3. To establish the effect of temperature on the UNIQUAC binary interaction 

parameters 

4. To compare the performance of different solvents towards their extractability for 
aromatics by using estimated parameters for UNIQUAC. 

5. To estimate the effects of aromatics content of feed, temperature, and solvent to 
feed ratio on selectivity and solvent power. 

6. To categorize the different solvents according to there performance for aromatic 
extraction. 
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Chapter 2 

LITERATURE REVIEW 


2.1: Solvent Selection Strategies 


As compared to distillation, economics of liquid-liquid extraction is not well 
stabilized. The economical operation of liquid-liquid extraction largely depends on 
selection of solvent. A desirable solvent must have a high selectivity (relative 
separation between two components) for one hydrocarbon relative to the other. The 
selectivity of the solvent is determined by the activity coefficients of the hydrocarbons 
when dissolved in the solvent. 

2.1.1 Limiting activity coefficients 

Infinite dilution activity coefficients play an important role in screening of 
potential solvents for liquid-liquid extraction. The solvents for aromatic recovery by 
either extractive distillation or extraction can be neatly rated in terms of the limiting 
activity coefficients for hexane and benzene as shown by ( Deal and Derr, 1964 [19] ). 
Limiting activity coefficients for hexane and benzene have been used by them to 
interrelate the capacity and selectivity for hydrocarbons for a broad series of solvents at 
25 to 100 °C. Selectivity has been found to decrease with increasing temperature for 
almost all solvents. ( Deal and Derr, 1964 [19] ) have reported Sulfolane, Dimethyl- 
Sulfoxide and Dimethyl-Formamide as potential solvents for aromatic extraction. 

A correlation for calculation of activity coefficients for different hydrocarbons 
when dissolved in excess of solvent has been given by ( Weimer and Prausnitz, 1965 
[82] ). This correlation can be used for screening of the solvents when no experimental 
LLE data are available, as it requires only pure component properties. But this 
correlation has a limitation because of its dependency on regular solution theory. 
Because regular solution theory can not predict activity coefficients less than one, the 
use of this correlation is limited. 

2.1.2 Industrially used aromatic extraction processes 

To compare the solvent performance for aromatic extraction, it will be 
beneficial to consider the existing aromatic extraction processes, apart from searching 
for new solvents. Extraction properties of Dimethyl-Sulfoxide (DMSO) as a solvent for 
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aromatic extraction have been discussed by Choffe et al., 1966. Extraction at ambient 
temperature using DMSO has been suggested operationally practical by them. This 
approach contrast against high temperature extraction processes (TETRA using 

Tetraethylene glycol, Udex using Diethyleneglycol water mixture etc.) used industrially 

% 

(Hydrocarbon processing, 1970). Extraction properties of N-Formyl-Morpholine (NFM) 
used by FORMEX process has been discussed by Cinelli et al., 1972. Selectivity and 
capacity of NFM based on infinite dilution activity coefficients have been compared with 
that of NMP and Sulfolane. The selectivity of NFM have been reported better then that of 
NMP. Solvent capacity of NFM has been reported higher than that of Sulfolane, but the 
selectivity has been reported poorer than that of Sulfolane. Effect of water on solvent 
performance has also been discussed and it was found that presence of water increases 
selectivity, while it decreases the capacity. Extractability of Sulfolane and DMSO along 
with other sulfur group solvents (2-4-Dimethyl Sulfolane, Vinyl Sulfone, Tetramethylene 
Sulfoxide, Glycol Sulfite etc.) has been compared based on limiting activity coefficients 
by Rawat et al., 1976. Solvents having sulfone and sulfoxide groups were found to be 
better as compared to other sulfure group solvents. A comparison of existing aromatic 
extraction processes has been made (Bailes et al., 1976). Sulfolane has been found to 
posses largest selectivity value and NMP the largest capacity value. Extraction properties 
of Sulfolane used by Shell Process have been discussed in (Hydrocarbon Processing, 
1982). Table 2.1 taken from Lo and Baird, 1980 gives an overview of existing ^omatic 
extraction processes. This Table also gives operating conditions along with the solvents 
employed. Sulfolane used as a solvent by Shell process gives good selectivity and a 
favourable solvent capacity and thus leads to low solvent to feed ratio required. Initially 
the glycol used by Udex process was diethylene glycol with water, but later it was 
replaced by tetraethylene glycol and water mixture. Tetraethylene glycol and water 
mixtures are claimed to increase solvent capacity by a factor of 4 and also require no 
antifoaming agent. Tetraethylene glycol is also used in TETRA process of Union carbide, 
but without mixing with water. A two stage extraction process is used in the TETRA 
process. DMSO process of Institut Francais du Petrole (IFP) uses dimethyl-sulfoxide 
along with 2% water as solvent. High density and low freezing point of DMSO makes it 
economical to extract aromatics at ambient temperature. N-methyl-2-pyrrolidone of 
Arosolvan process is generally blended either with ethylene-glycol (40-50 wt %) or water 
(12-20 wt %) to increase the selectivity and to decrease the boiling point of the solvent. 
However addition of water 
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or ethylene- glycol slightly reduces the solvent capacity. The quantity of co-solvent 
required depends on the aromatics content of the feed. 


Table 2.1. Processes used for extraction of aromatics 


Process 

Solvent 

Operating 

conditions 

Comments 

Shell Process 

Universal Oil 

Products 

Sulfolane 

120 °C 

High Selectivity and 

capacity of sulfolane leads 

to low solvent to feed ratio 

Udex Process, 

Universal Oil 

Products 

Diethylene- 

glycol/water 

mixture 



o 

o 

O 

Tetraethylene glycol and 

water mixtures are claimed 

to increase capacity by a 

factor of 4. A two step 

distillation to recover BTX 

Union Carbide 

Corporation 

Tetraethylene- 

-glycol 

(TETRA) 

o 

o 

o 

O 

The extract leaving the 

primary extractor is 

essentially free of feed 

aliphatic, and no further 

purification is necessary 

Institut Francis 

du Petrole 

Dimethyl- 

sulfoxide 

(DMSO) 

Ambient 

Solvent has low freezing 

point and is nontoxic and 

noncorrosive 

Arosolvan 

Process, Lurgi 

N-Methyl-2- 

Pyrrolidone 

(NMP) 

NMP-ethylene - 

glycol (40-50 wt%), 

60 °C 

Either ethylene-glycol or 

water is used as cosolvent 

with NMP 

The quantity of co-solvent 
required depends on the 

aromatics content of the 

feed. 

NMP-water 

water( 12-20 wt%) 

35 °C 

Formex Proess, 

SNAM Progetti 

[30] 

N-Formyl- 

Morpholine 

(NFM) 

40 °C 

Low corrosion allows use of 

carbon steel equipment 
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2.1.3 Molecular Design of Solvents 

On the basis of the UNIFAC group contribution, it was proposed to synthesize 
molecular structures with specific solvent properties for the separation of aromatic and 
paraffinic hydrocarbons by ( Gani and Brignole, 1983 [27] ). The method of molecular 
design of solvents has been well formulated in terms of combination rules and target 
properties ( Gani et. ah, 1991 [26] ). It has been observed that large number of 
molecules can be generated by considering small number of groups. The size of 
combinatorial problem posed by the molecular synthesis procedure is reduced by 
group selection and by imposing physical and molecular constraints at different 
stages by ( Pretel et. ah, 1994 [58] ). Limitations of UNIFAC arise due to non- 
availability of binary interaction parameters for groups. To circumvent this problem 
(Meiiiai and Newsham, 1995 [42]) has proposed method of calculating group binary 
interaction parameters, necessary for molecular design, using molecular graphics 
system. 

2.2: Parameter Estimation 

For evaluation of solvent performance for liquid-liquid extraction, the required 
parameters (selectivity, solvent power etc.) are calculated through prediction of liquid- 
liquid equilibria. For prediction a suitable thermodynamic model which can represent 
liquid-liquid equilibria should be used ( Bailes et. ah, 1976 [7]). The excess free energy 
models such as NRTL the local composition model ( Renon & Prausnitz, 1968 [67] ) 
and the quasi chemical theory models UNIQUAC ( Abrams & Prausnitz, 1975 [1] ) and 
UNIFAC ( Fredenslund et ah, 1975 [25] ) can represent liquid-liquid equilibria fairly 
well. These models require proper binary interaction parameters, which are generally 
obtained by regression of experimental data. Prediction of liquid-liquid equilibria is 
very much sensitive to the choice of parameters as compared to the vapour-liquid 
equilibria ( Abrams and Prausnitz, 1975 [1] ). For estimation of binary interaction 
parameters, binary LLE data for partially miscible systems and binary VLE data for 
completely miscible systems can be used. But for prediction of ternary liquid-liquid 
equilibria, which is normally used in liquid-liquid extraction calculations, the 
parameters obtained by using this approach have been found to be poor ( Varhegyi and 
Eon, 1977 [78]). This is mainly due to the fact that reduction of binary VLE data for 
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completely miscible systems cannot yield imique results unless the accuracy of the data 
is extremely high, much higher than usually encountered in typical experimental 
studies. Within the experimental error, many sets of binary parameters can equally well 
reproduce the binary VLE data (Abrams and Prausnitz, 1975). Furthermore if the 
binary interaction parameters thus obtained are used for multicomponent LLE 
prediction then the errors may be too large to justify (Anderson and Prausnitz, 1978). 
To avoid these discrepancies Varhegyi & Eon, 1977 have shown that ternary LLE data 
should be employed for regression of binary interaction parameters, in this way effect 
of VLE can be avoided. They have outlined the procedure to regress binary interaction 
parameters for NRTL model using ternary LLE data. Minimization of square of 
weighted errors between experimental and predicted compositions has been used to get 
the best fit. Varhegyi & Eon make no attempt to discern the physical meaning of the 
calculated parameters. The parameters have been treated only as correlation 
coefficients. To avoid the problem of getting different parameters for same pair, (Chen 
et al., 2000) have suggested to keep the parameters for a particular pair to be same for 
all systems. 

A basic assumption of the least-square minimization is that the values of the 
independent variables are known without error. This assumption fails to account for the 
statistical properties of the independent variables, arbitrarily assigning them a zero 
standard deviation; therefore it does not utilize all available information in estimating 
the parameters. A more suitable method, which takes into account the errors in 
independent variables, will be the one based on maximum likelihood principle as 
suggested by Anderson et al., 1978. They have used maximum likelihood principle for 
the determination of model parameters from experimental data when all the measured 
variables are subject to error. Method has been applied for reduction of binary vapour- 
liquid equilibrium data. Sorensen and Arlt, 1979, have used maximum likelihood 
approach for regression of binary interaction parameters using ternary LLE data. 
Almost all of the maximum likelihood approaches currently used to regress the 
parameters of thermodynamic models fix the variances, converting the problem into a 
traditional weighted least square minimization. The variances are generally fixed on the 
basis of experimental error estimates (i.e. standard deviations). However such 
approaches lead to residual variances (between measured and calculated composition 
values) that are inconsistent with the fixed variances and thus, do not produce optimum 
parameters for prediction purposes. ( Vasquez & Whitting, 2000 ) 
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have recently suggested a method, referred as IVEM (Inside Variance Estimation 
Method) which involves the re-computation of variance for each iteration of the 
optimization procedure, automatically re-weighting the objective function. This results 
in better parameter estimation and to a direct measure of uncertainty in the model 
prediction. In a later article (Vasquez & Whitting, 2000 [79] ) has shown the effect of 
regression approaches in the estimation of nonlinear model parameters on the predicted 
process performance. 
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Chapter 3 

Solvent selection strategies 

3.1 Important characteristics of solvents for aromatic extraction 

An efficient solvent for aromatic extraction should have sufficient solubility for 
aromatics, while very low solubility for non-aromatics. The solubility is mainly 
influenced by activity coefficients. Apart from solubility, the following properties of 
solvents, which can be characterized into primary and secondary properties, should be 
considered to evaluate the performance of solvents for extraction of aromatics. 

Primary properties 

> High Selectivity for aromatics 

> High Capacity for aromatics 

> Low Solvent Loss 
Secondary properties 

> High Density difference with Feed 

> High Boiling Point 

> Thermal Stability 

> Sufficiently High Freezing Point 

> Non Toxicity 

>■ Environmentally acceptable 

> Non Corrosiveness 

> Low Price 

Among these the first two, i.e. solvent selectivity and capacity are the 
two primary factors, which alter the efficiency of extraction process. For solvent 
selection these two parameters are considered primarily, with more importance given to 
selectivity. However, the best solvent is not necessarily that with highest selectivity 
value, since high selectivity often means low solvent capacity or solvent power for the 
extracted component. Solvent capacity has an even greater influence than selectivity on 
the selection of optimum extraction process because the capacity determines the flow 
rate of the circulating solvent and consequently the size of extractor. These two primary 
factors are defined as follows: 
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3.1.1 Solvent Selectivity 

Selectivity is defined as the ratio of the distribution coefficients of the one 
aromatic component to that of one non-aromatic component in the system of solvent - 
aromatics - non-aromatics. 

Hence, 

, . . Distribution coefficient for an eromatic 

selectivity = ^ 

Distribtution coefficient for a non-aromatic 


or, in terms of notations (see nomenclature) 



x/ x" 


(3.1) 


where subscripts ‘i’ and ‘j’ refer to the aromatic and non-aromatic components and 
superscripts I and II denote the extract and raffinate phases, respectively. Following the 
thermodynamic criteria for the phase equilibria, described later, it can be shown that 
(Mukhopadhyay, 1979 [46] ): 



rirj 

/y; 


(3.2) 


3.1.2 Solvent Capacity 

Solvent capacity or the solvent power is defined as the distribution coefficient of the 


solute (aromatics). 



(3.3) 


following the thermodynamic criteria for phase equilibrium, it can be written as 


5 


P 




(3.4) 


It is usually observed that the solvents having high selectivity will have low 
solvent capacity and vice versa. Therefore a compromise between these two parameters 
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is necessary, for which a solvent selection criterion which will incorporate both these 
parameters is required. A simple criteria, the product of Selectivity and Solvent 
Capacity, called as performance index (PI) and defined as: 

PI = 5*5,, (0.1) 

was suggested by Brignole et al., 1986. However more terms such as Solvent Loss, 
Toxicity, etc. can be incorporated in the definition of PI depending on the particular 
needs. 

3.2 Search for new solvents 

Although most of the solvents commercially used possess good selectivity, their 
capacities are low. Therefore a high solvent flow rate is needed to achieve the desired 
recovery of aromatics. Because the extraction processes are found to be very expensive, 
it will be economical to search for new solvents, which can have favourable capacity 
apart from good selectivity. Good methodologies have been developed to search for 
new solvents for extraction as well as extractive distillation. The methodologies known 
as molecular design of solvents are widely used for development of solvents for 
separation processes. 

3.2.1 Molecular Design of Solvents 

Molecular design of solvents, originally developed to select entrainer for azeotropic 
distillation, has been extended for extraction solvent selection by Gani & Brignole, 
1983. In this methodology the solvent molecules are designed from their constituent 
groups, so that the solvent will have the properties to achieve the desired separation. 
The solvent molecules are designed, keeping in view the desired properties such as 
selectivity, solvent power etc. These methodologies, primarily discussed by Gani & 
Brignole, 1983, Brignole et al. 1986, Gani et al, 1991, Pretel et al., 1994 and 
Constantinou et al, 1996 are based on the group contribution approach UNIT AC. In the 
UNIT AC method, the molecules are decomposed into groups and the properties of 
components are then calculated based on the group properties. A reverse technique is 
used for molecular design of solvents-here the molecules are generated with groups, 
which can have the desired properties. The overall procedure for molecular design of 
solvents can be broken into following steps: 

(a) Classification of groups 

' The groups, which are to be considered for molecular design of solvent, are broadly 
categorized as terminal or intermediate groups ( Brignole et. al., 1986). Groups, 
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having single attachment for combination, are called terminal and having multiple 
attachment are called intermediate groups. These groups are further classified 
according to classes and categories ( Gani et. al, 1991 [28] ). The class of the group 
signifies the number of free attachments available to the group. The category of a class 
of groups signifies the degree of restrictions they have with respect to their linkage with 
other groups. 

(b) Selection of groups and target properties 

The target properties are selected according to the separation process for which, the 
solvent is to be designed. For extraction solvent design, the target properties may be, 
the selectivity, solvent power, solvent loss, or any combination of these. After selection 
of target properties, a sensitivity analysis of these properties with respect to all the 
groups is conducted, to pre-select not only a sub-set of groups but also, to pre-select 
certain type of compounds (aromatics, alcohols, ketones etc.). This analysis attempts to 
determine the most sensitive groups (called functional groups). The sensitivity analysis 
indicates the groups which are most likely (or least likely) to be present in the best 
molecular structures. The interchange energy, (IE), defined by Magnussen et al.,1981 
[40] as an average of the binary interaction parameters, is used as a criteria for the 
selection of groups for a specific problem. A large positive value for the IE means 
solubility for the solute in the solvent. The interchange energy between groups i and j 

can be calculated from a summation of the group - interaction parameters a;j and as 


given by Magnussen, 1981 [40] 

" R ’ R 

a,! + a = [u,j - ^ (w, + )] I = (/E). I 


(3.6) 

(3.7) 


A small negative value of (IE) indicates good affinity between the solute group and the 
solvent. The difference in the values of IE for two groups in a solvent gives an 
indication of selectivity. If the differences are small and the values are negative, the 
ternary system will exhibit complete miscibility and separation with the solvent will be 
infeasible. The prediction of solvent properties requires the availability of UNIFAC 
binary interaction parameters. Therefore, the first criterion for the selection of groups is 
the availability of binary parameters for the synthesis groups with those of the feed 
mixture compounds. 


14 



On the basis of above the groups that should give good solvents for aromatic 
extraction have been identified. These groups are given in Table 3.1(a) along with 
references. 


Table 3.1(a) Groups identified for aromatic extraction solvents 


S.N. 

Group 

Reference 

1 

CH^O , 

Menial &Newsham, 1999 [43] 

2 

CH,0, 

Gani & Brignole, 1983 [27] 

3 

CHfiN 

Brignole et. al., 1986 [9] 

3 

CHNO^, 


4 

CH^NO^, 

Gani & Brignole, 1983 [27] 

5 

CH^CO 

Brignole et. al.,1986 [9] 

6 

ACNH^, 


7 

CHfiOO, 

Gani & Brignole, 1983 [27] 

8 

ACCH, 



(c) Formation of feasible chemical compounds (combination rules) 

The generation of structures from the groups is subject to certain restrictions. These 
restrictions are imposed to prevent formation of unstable compounds or the appearance 
of new functional groups, like acetals, for which the prediction of the properties will be 
uncertain. The molecular structures to be formed must satisfy the primary condition 
that the final molecular structure cannot have any free attachments Groups of certain 
categories have restrictions for their joining with other groups; these restrictions must 
be fulfilled. These restrictions, commonly referred to as combination rules, are well 
developed and tabulated. 

(d) Evaluation of solvent properties and screening of compounds 

The primary properties of a solvent for liquid-liquid extraction are selectivity and 
solvent capacity. The secondary properties are boiling point, density and viscosity etc. 
Once the chemically feasible molecular structures are formed, they can be tested for the 
constraints on the target properties. At this stage therefore the property prediction 
strategies are needed. The group contribution method UNIFAC is useful at this stage. 
However the binary interaction parameters for all the constituent groups of the 
generated solvents with the groups in feed should be available. Therefore the 
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availability of binary interaction parameters of all groups which are to be considered is 
of vital importance. But these parameters are still not available for large number of 
groups, which hinders the efficient use of molecular design of solvents. 

A similar approach of molecular design of solvents developed initially by 
Meniai & Newsham 1992 [41], uses molecular graphics system for calculation of 
binary interaction energy for groups. A crude method of calculating binary interaction 
energy’s, based on molecular graphics system coupled with Monte Carlo simulations 
has been proposed by Meniai & Newsham, 1995 [42]. The method has been tested for 
solvent selection for two industrial extraction processes Meniai & Newsham ,1999 
[43]. Promising results have been obtained by this method. 

3.3 Screening of Solvents based on Infinite Dilution Activity 
Coefficients 

The activity coefficient at infinite dilution (limiting activity coefficient) represents 
an important property, which is used in particular for the selection of selective solvents 
and for the reliable design of separation processes. The study of infinitely dilute 
solutions deserves considerable interest from both practical and theoretical point of 
view. In fact if we consider a mixture in which an important component is present at a 
very low concentration, the knowledge of the behavior of the dilute solution is essential 
for evaluating the recovery and separability of the components. The activity coefficient 
at infinite dilution characterizes the behavior of the solute molecule surrounded by the 
solvent: this indicates a maximum non-ideality whereas at finite concentrations the 
non-ideality is “diluted”, that is the activity coefficients assume values closer to unity. 
Therefore to be able to estimate the selectivity of a solvent, it is often sufficient to 
know the separation factor at infinite dilution. The activity coefficients at infinite 
dilution also provides a measure of separability of different solvents. 

At infinite dilution, the parameters defined earlier for selection of solvents takes the 
form shown below. 

3.3.1: Selectivity at infinite dilution of solutes 



At infinite dilution the ratio { y 7 ^ in nonpolar phase, becomes nearly 
equal to unity, as recommended by Deal & Derr, 1964 [19]. Therefore the selectivity at 
infinite dilution of solute approximately becomes 
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(3.9) 


3.3.2: Solvent capacity at infinite dilution of solute 

The solvent capacity at infinite dilution of solute can be approximated by the 
relationship 



(3.10) 


This approximate relation holds for almost al systems, involving aromatic and paraffin 
hydrocarbons, as proposed by Deal & Derr, 1964 [19], 

3.3.3: Performance Index at infinite dilution of solute 


On the basis of above definitions of solvent selectivity and capacity, the performance 
index defined in equation (3.5) for comparison of solvent performance can be 
approximated as follows. 



(3.11) 


These infinite dilution activity coefficients are generally obtained 
experimentally. The differential ebulliometry and gas liquid chromatography are the 
two most widely used techniques for experimental measurements of infinite dilution 
activity coefficients, with later gaining more importance now a days. If experimental 
value is not available for any promising solvent, it is possible to calculate these by pure 
component data only. The correlation’s based on regular solution theory, such as one 
developed by Weimer & Prausnitz, 1965 [82], are useful for this purpose. The 
correlation developed by Weimer & Prausnitz, 1965 [82] is applicable for calculation 
of infinite dilution activity coefficients, for paraffin’s, naphthenes and aromatic 
hydrocarbons in a large variety of polar solvents. Other generalized correlation’s such 
as MOSCED (modified separation of cohesive energy density) gives good estimate of 
activity coefficient at infinite dilution ( Reid et. ah, 1988 [66] ). 

3.4: List of potential solvents and screening 

Based on the above strategies and on the basis of industrially used solvents for 
aromatic extraction processes a list of potential solvents has been made. This list also 
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includes the solvents suggested by the method of molecular design of solvents reported 
by different authors. 

These solvents have been screened based on their performance index at infinite 
dilution (PI“). For calculation of selectivity, solvent capacity and performance index 
calculation the experimental values of activity coefficients at infinite dilution (y“) are 
used. The effect of temperature on selectivity and solvent capacity can also be looked at 
based on the variation of infinite dilution activity coefficients with temperature. The 
effect of temperature on selectivity, solvent capacity and performance index are given in 
Figure 3.1, Figure 3.2 and Figure 3.3 respectiviely. The references of the y“ values used 
for calculation of selectivity and solvent capacity are reported in Table 3.1(b) 


Table 3.1(b) Performance Index for Solvents along with References 


S.N. 

Solvent 

S” 

SP“ 

(pr) 

Ref. for 

00 

Y 

1 

N-Methyl-2Pyrrolidone 

(NMP) 

ll59 

0.96 

12.11 

[81] 

2 

Sulfolane 

(TMS) 

19.76 

0.42 

8.31 

[44] 

f ‘I*?!* 

3 

N-Formyl-Morpholine 

(NFM) 

16.01 

0.45 

7.80 

[81] 

He-iJCich, 

4 

Dimethyl-Sulfoxide 

(DMSO) 

20.36 

0.30 

6.17 

[75] 

Ti«^s , n 8 6’ 

5 

Dimthy 1-F ormamide 

(DMF) 

12.12 

0.44 

5.30 

[76] 

6 

T etraethy lene-Glycol 

(TEG) 

13.80 

0.38 

1 

5.28 

. 

[75] 

7 

T riethy lene-Glycol 

(TREG) 

17.50 

0.22 

3.95 

[75] 

■ 

Propylene-Carbonate 

(PC) 

12.19 

0.30 

3.66 

[77] 

9 

Diethylene-Glycol 

(DEG) 

13.00 

0.12 

1.70 

[28] 

me 

10 

N -Methy 1-F ormamide 

(NMF) 

f 2- 

0 • l-J 

1.64 

[30] 
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Fig.3.1 : Variation of Selectivity in excess solvents with temperature 
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Fig. 3.2: Variation of solvent capacity in the excess solvent with temperature 
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Fig.3.3: Variation of Performance Index in excess solvents with teperature 
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Chapter 4 

PARAMETER ESTIMATION 

The evaluation of solvent performance for analysis of aromatic extraction 
systems requires huge amount of experimental data. But experimental data on multi- 
component, (i.e. involving more than three components) liquid-liquid equilibria (LLE) 
of aromatic extraction systems are rather scarce. Therefore reliable prediction of 
multicomponent data from the knowledge of limited experimental data information is 
of vital importance to accurately evaluate the performance of solvents. Prediction of 
multi-component LLE and parameters for solvent performance evaluation requires the 
knowledge of activity coefficients in terms of compositions (Mukhopadhyay and 
Dongaonkar, 1983). Several models, such as UNIQUAC, UNIT AC, NRTL have been 
proposed to represent the composition dependence of the activity coefficients in 
mixtures (Walas, 1985). Among these three models the UNIQUAC model having two 
adjustable parameters gives almost similar predictions as compared to NRTL model 
having tlnee adjustable parameters (Prausnitz et. al, 1980). All these models require 
binary interaction parameters. Liquid-liquid equilibria are found to be so sensitive to 
the value of tlrese binary parameters that a small uncertainty in parameters can produce 
drastic errors in prediction of compositions (Gomis et. al., 1997). These model 
parameters are estimated by regression of experimental liquid-liquid equilibrium data. 
Small errors in experimental data can produce large uncertainty in regressed 
parameters, which in turn may lead to unreliable prediction (Vasquez & Whitting, 
2000). Therefore the regression technique chosen should be capable of handling 
random errors in the experimental data. (Xin et al., 2000) has shown the effect of 
parameter regression approach on the prediction of liquid-liquid equilibria. 

Our purpose here is to get thermodynamic model parameters, which can give 
accurate prediction of activity coefficients and compositions in nonideal liquid 
mixtures. Therefore for parameter estimation we should keep in view the minimization 
of errors between experimental and predicted compositions (or activities). Hence the 
regression of parameters is a minimization problem, from a mathematical point of view. 
The two main approaches to define objective function for minimization are least 
square minimization and maximum likelihood principle. Out of these two general 
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techniques, maximum likelihood approach takes into account the possible random 
errors in experiments (Prausnitz et. al., 1980 [56] ). However the least squares approach 
needs less computational effort. 

Regression of binary interaction parameters needs predictions of composition 
variables for objective function calculation. Therefore the regression of liquid-liquid 
equilibria should start with solution of phase equilibria equations. The necessary 
condition for phase equilibria 

LLE: y\ xj = rl'x'/ , i =1,2 c (4.1) 

along with summation equations, 

c c 

Summation: (4.2) 

(=1 1=1 

are solved at each iteration. The equation (4.1) is valid only when, same standard state 
fiigacity is used in both the liquid phases. An acceptable assumption is that pressure, 
has no effect on the liquid-liquid equilibria, so the activity coefficients are functions 
only of temperature and compositions. 

After prediction of phase equilibria using a suitable initial guess of parameters 6, the 
objective function is calculated next. 

4.1 Least square minimization 

Representing the parameters to be regressed by 9 vector, the most simple 
objective Imnction using the least square technique, has the form 

m c 2 

1 (0) = min X (4 > 6) " 4 4 (4 ^ J (4.3) 

In this objective function the summation of squared differences in activities of 
components in two phases is minimized. Many variations of the general objective 
function represented by equation (4.3) are possible, as given by ( Vasquez and Whiting, 
2000 [79] ). Using the same approach, least square minimization of activity differences 
may take the form 

m c _2 

(0) = min 2] [in 4 (x^ , 0) - In 4 (xjl , 0) J (4.4) 

® k=\ i=l 
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This objective function is a relative objective function, in which the ratio of activities is 
used. The objective function of equation (4.3) is an absolute objective function. These 
activity based objective functions have the common problem of not yielding good 
agreement with experimental results in general ( Anderson et. al., 1978 [2] ). However, 
they have very good convergence properties. 

In the estimation of binary interaction parameters, our aim is to get the set of 
parameters for thermodynamic model, which can give good prediction of compositions 
of a multi-component mixture in both phases. Therefore more appropriate objective 
function should be the one which minimizes the errors between experimental and 
predicted compositions. The most general objective function of this type, with a least 
square minimization technique is 


'^j(9) = minj ('t-S) 

k = [ j=\ /=1 

This type of objective function requires prediction of phase compositions, at each 
iteration of regression. For phase composition prediction and binodal curve prediction, 
the phase equilibria equations (4.1) along with summation equations (4.2) are solved 
simultaneously. And this prediction of phase equilibria is found to be computationally 
intensive and needs a qualified initial guess of parameters. Because of their good 
convergence, the activity based objective functions can be used at initial stage of 
regression procedure. The parameters obtained by minimization of objective functions 
or initial guess for the objective function of equation Fjj. 


Alternatively any other initial guess of parameters which can predict both phases can be 
used. 

A basic assumption of least square method is that it discriminates between 
independent and dependent variables during regression. This method assumes that the 
independent variable (measured variable) is known without errors and only the 
dependent variables are susceptible to error. But this is not the case generally for the 
parameter estimation from liquid-liquid equilibrium data. Because the measured 
variables, i.e. the experimental compositions are always found to contain some errors 
(Prausnitz et. al., 1980 [57]), least square method cannot give the best estimate of 
parameters. This method does not account for statistical properties (standard 
deviations) of the independent variable. Since almost all models used for liquid-liquid 
equilibria are nonlinear in their adjustable parameters, the best estimates of these 
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parameters can be obtained from a formalized method, which properly treats the 
statistical behaviour of the errors associated with all experimental observations. 

The methods based on maximum likelihood principle attempts properly to 
account for all measurement errors, as discussed by Anderson et al, 1978. 

4.2 Maximum likelihood estimation 

It is usually assumed that there are two types of measurement errors; systematic 
and random. The formers are due to inherent bias in the measurement procedure 
resulting in a consistent deviation of the observable from its true value. Usually, an 
experimenter’s skill and experience provide the only means of consistently detecting 
and avoiding the systematic errors. Random or statistical errors are assumed to be the 
result of a large number of small disturbances. Such errors have been found to be 
distributed according to simple laws that make it possible to treat them with statistical 
methods. A major consequence of these random errors is the corresponding presence of 
errors or uncertainties in the estimated parameters. Because of these errors in data, and 
also because of simplifications and assumptions in the model, it is not possible for a 
model to represent the experimental data exactly. However, a method of parameter 
estimation, which correctly utilizes all the pertinent information available will give a 
best fit of the model to the data and minimize the parameter uncertainty. The principle 
of maximum likelihood provides such estimation. 

4.2.1 Maximum likelihood principle 

In the maximum-likelihood analysis, it is assumed that all measured data are 
subject to random errors. If the experiments were replicated, the average value for each 
replicated experimental point would approach some true value. Usually the distribution 
of a measured variable about its true value is approximated by the normal distribution, 
characterized by an associated variance. If there is any coupling between the 
measurement method, then there are also associated covariances between the measured 
variables. These variances and covariances must be known or estimated, although 
covariances are almost always assumed to be negligible (Anderson et al., 1978). It is 
customary to assume that the random errors in the different experiments are 
uncorrelated. 

For each experiment, the true values of the measured variables are related by 
one or more constraints, those that account for material balance and phase equilibria. 
Because the number of data points exceeds the number of parameters to be estimated. 


25 



all constraints are not exactly satisfied for all experimental measurements. Optimum 
parameters and the true values corresponding to the experimental measurements must 
be found by satisfaction of an appropriate statistical criterion. 

If this criterion is based on the maximum-likelihood principle, it leads to those 
parameter values that make the experimental observations appear most likely when 
taken as a whole. The likelihood function is defined as joint probability of the observed 
values of the variables, model parameters, and error variances. The best estimate of the 
model parameters and of the true values of the measured variables are those, that 
maximize this likelihood function with a normal distribution assumed for the 
experimental errors. 

4.2.2: Maximum likelihood estimation 

The model for maximum likelihood parameter estimation takes the form of an 
implicit algebraic system of equations 

f(0,z) = 0 (4.6) 

where 0 is a vector of unknown parameters, z is a vector of n experimentally measured 
variables, and f represents the system of algebraic functions. 

The measurements are related to the true values through 

Za = m (4.7) 

where ‘ m ’ is number of experimental data ( tie lines), and is vector of unknown true 
values of the experimentally measured variables, , at the A:* tie line, and is vector 
of errors. 

The likelihood function (L) is defined as the probability of the observed errors 
in all data points occurring given a set of statistical parameters {^/), i.e. the standard 
deviation for the distribution used to represent those errors. As it is generally assumed 
and is almost valid (Prausnitz et. al, 1980 [57]), the measurement errors from different 
data points are considered as uncorrelated or independent. Then, the joint probability 
density can be expressed as a product of the individual densities. Hence the likelihood 
function should be (Esposito & Floudas, 1998 [21]) 

m 

Li^,¥) = ^P{^ki^)\¥k) ( 4 - 8 ) 

where p represents the probability density for the individual tie lines. 

Further if we assume that the errors for n measured variables follow a normal 
distribution with zero mean and covariance matrix V, we have 
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Here « = 2 for a binary system, n = 4 for ternary system, and for a multi-component 
system n = 2*(c-l), where c is the total number of components. 

Hence, from equations (4.8) and (4.9) the logarithm of likelihood function L can 
be written as 


nnj 1 1 

In i(9, (,) = - - -In 2;r --y; In I V. I - - 2; e, V;‘e, (0.2) 

^ ^ k=\ ^ i=l 

The maximum likelihood estimator seeks to determine values for 0 and 
^ which maximize L. In the estimation of binary interaction parameters for 
thermodynamic models it is always assumed that the errors in each experiment are 
known, equal, and independent, that is, the covariance matrix V, is diagonal with 
elements v,.; , then the maximization of L from equation (0.2) is equivalent to 

m n 

F 2 , (0, y/) = min 2 ,.v:’ (0.3) 

A:=l /=! 

where e ,^ , is the i”' component of the vector e^. , that is the error associated with the i'^' 
variable in the k"' experiment. 

This is the most widely used objective function for parameter estimation of nonlinear 
thennodynamic models by maximum likelihood function (Sorensen & Arlt, 1979). In 
this approach an priori estimate of variances are made based on experimental error, thus 
converting it to weighted least squares minimization. The variance-covariance matrix is 
fixed based on experimental precision. But as shown (Vasquez & Whiting, 2000), this 
method of using an priori estimate of variances is internally inconsistent because this 
estimate is not (in general) equal to the final values obtained from the residuals at the 
conclusion of the regression procedure. To avoid this inconsistency Vasquez and 
Whiting, 2000, suggested a method known as Inside Variance Estimation Method 
(IVEM). Although this method uses the objective function given by equation (4.10), 
V;; is estimated from each iteration in the regression procedure. Thus, at the end of the 
optimization process, the statistical properties will be consistent with V^ , guaranteeing 
that the most likely values of the 0 are found. In this method also it is assumed that 
errors are independent with each other, 


27 



so that the variance-covariance matrix is diagonal, with zero covariance terms. 
Therefore the maximization of L of equation (4.10) can be written as minimization of 

F 22 (0, y/) = min(- In L) = min | ^ In 2;r + — ln(J^ v,.^) + - ^ V;'e^ | (4.12) 

I 2 2 2 J 

The regression process starts with a initial guess of parameters, which is used 
for phase envelope generation, and the objective function of equation (4.12) is 
minimized with any appropriate optimization technique. The predicted compositions 
are then compared with experimental ones either by objective function value or using 
root mean squared deviation {rmsd) defined as, for a ternaiy system 


rmsd = 1 00 


m n 2 

TIE 


^ = 1 / = 1 ,/=:l 


6m 


1/2 


(4.13) 


here 6 in the denominator indicates the number of composition variables, 3 in each of 
the two phases. For a binary system it will be 4. 

The algorithm for estimation of binary interaction parameters for thermodynamic 
models using liquid-liquid equilibria can be written as follows. 

4.3: Algorithm for parameter estimation using LLE data 

The main algorithm can be divided into following steps: 

Step 1 : Initial guess of parameters 

A qualified initial guess of parameters, which is crucial for regression of LLE, 

can be obtained from different sources such as 

>• Aspen plus process simulator: Use parameter estimation run of Aspen plus to 
get UNIQUAC model parameters, using predictive UNIFAC for LLE 

> Activity difference approach: Minimize objective function of either equation 
(4.3) or (4.4), with any initial guess. 

> DECHEMA data bank: Interaction parameters for large number of systems 
have been reported in DECHEMA data bank (Sorensen & Arlt, 1980). Use 
common parameters as initial guess. 

> Reported LLE parameters: Use any reported values of parameters from any 
other source in the literature as initial guess. 

> Binary VLE parameters: If no liquid-liquid parameters are available, use 
binary VLE parameters for miscible pairs and binary LLE parameters for 
partially miscible pairs, as initial guess. 
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step 2: Predict phase equilibrium compositions (tie-lines) 

Use the initial guess of parameters obtained in step 1 to estimate activity 

coefficients using the thermodynamic models of concern 
(UNIQUAC/UNIFAC). Activity coefficients as a function of binary interaction 
parameters, composition and temperature are given in Appendix-A for 
UNIQUAC and UNIFAC models. 

Store experimental tie lines into matrix est_tielines (m, n) 
for all experimental tie-lines exp_tielines (m, n) do 

Compute the estimated tie line corresponding to the experimental one. 
Estimate the tie line, which passes through mid-point of experimental 
tie line. For this solve the simultaneous set of equations given below 
For binary systems: Solve 

Phase equilibria: x’y' = x'/yf i = 1,2 

along with 

2 2 

Summation: ^x/=l and =1 

(=1 /=) 

for 4 composition variables, 2 in each of the two phases . For 
calculation of activity coefficients use the relations given in Appendix- 

A. 


(4.14) 

(4.15) 


For ternary systems: Solve 
Phase equilibria: 

along with 


^ / / ^ II II 

hr. =h Yi 


1 = 1,3 


Summation: 

and 

Predicted tie line: 


E^/=l and^xf=l 

/=] /=! 


(x, -X,) 


(4.16) 


(4.17) 


(4.18) 


for 6 composition variables, 3 in each of the two phases. Equation (4.18) 
is equation of tie line, which passes through the mid-point of the 
experimental tie line as shown in the figure. Where z, and Zj are the 
compositions at the mid point of experimental tie lines as shown in 
Figure 4.1, and given by 
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end 


(x. 


, +x") 


and 


{x[ +x") 
2 


if 

above equations are solvable and give real solution of composition 
variables for all tie lines, then 

store estimated tie lines in matrix est_tielmes(m, n), go to step 3 

else 

change the initial guess of parameters, go to step 1 

end 

Step 3: Evaluate objective function 

Evaluate the errors between exp_tielines (m, n) and est_tielines (m, n). 

(0) = exp_tielines(A:, /)-est_tielines(A:, /) 

k = 1, 
i = l,....,n 

Take n = 2 for binary systems and « = 4 for ternary systems 

for all i = 1 ,....,« do 

compute the average error for the z'* component over all tie lines 

m 

(4.19) 

m 

and the diagonal elements of variance in the matrix V given by 


V, =M 

" (m-1) 

assign zero values to the off diagonal elements of V 


(4.20) 


end 


Evaluate the objective function , given by equation (4.12). 


Step 4: Minimize objective function 

Minimize the objective function according to the optimization technique used. 

if 

objective function has reached a minimum then store parameters (0) as 
the optimum parameters. 
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else 

change 0 according to optimization procedure used, until minimum of 
objective function is reached. 

end 

Calculate the root mean squared deviation defined by equation (4.13) 

Although the algorithm is written for inside variance estimation method (IVEM) 
based on maximum likelihood function, it can be used for least squares minimization 
technique also, with minor changes in step 3. For least squares technique the objective 
function is calculated directly after step 2 without calculating other statistical terms. 
The same algorithm can be used for parameter estimation of UNIQUAC as well as 
UNIFAC. 

For minimization any suitable optimization technique can be used. The 
technique, used here is the simplex method, suggested by Nelder and Mead, 1964 [54]. 
Because this method does not require function derivative calculations, the round off 
errors will be low for this method (Anderson et. al., 1978 [2]). 
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Chapter 5 

Parameter Estimation Results 

The solvent performance for different solvents is evaluated on the basis of 
parameters such as selectivity, solvent capacity, solvent loss etc. as described in 
previous chapters. These parameters are defined in terms of activity coefficients, which 
in turn are calculated through thermodynamic models. The activity coefficient models, 
UNIQUAC and UNIT AC are widely accepted for activity coefficient prediction. These 
models are valuable only if the proper binary interaction parameters are available. 

The feed, from which the aromatics are extracted, contains large number of 
compounds and several families (Paraffin, Isoparaffin, Olefin, Naphthene, Aromatic). 
The components present in the feed covers a large range of molecular weights. 
I’herefore it is impossible to model exactly the aromatic extraction system. Furthermore 
the quantity of total aromatic and non-aromatic compounds present in the feed is also 
found to be different for different sources of naphtha. For evaluation of performance of 
solvents, the feed is generally assumed to be composed of key components. The 
number of key components used is also tentative and depends largely on the availability 
of experimental data. For the comparison of extractability of solvents, the feed is 
modeled as having one key aromatic component and one non-aromatic key component. 
More components can be added, based on the availability of experimental LLE data. 
The binary interaction parameters for the feed and solvent systems are than estimated. 
The systems used for parameter estimation here are given in Table 5.1(a). The 
multicomponent systems used for comparison of predictions using estimated 
parameters are listed in Table 5.1(b). 

5.1. Separate Estimation of parameters 

For the estimation of binary interaction parameters the ternary data alone can be 
used as suggested by (Verhegyi and Eon, 1977 [78]) and implemented in DECHEMA 
data book (Sorensen and Arlt, 1979 [71]). However, the parameters obtained by this 
technique are found to be specific to that particular ternary system from which they 
were regressed. For example at the same temperature the interaction parameters for 
Hexane-Benzene pair have been found to be different in the system Hexane-Benzene- 
Sulfolane and Hexane-Benzene-Dimethylsulfoxide, as reported in DEHCEMA. Similar 
behavior has been observed for all other systems also. 
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Table 5.1(a) Ternary Systems used for parameter estimation. 


System 

Number 

Ternary system 

Temperature (K) 

1 

Hexane - Benzene- TMS 

298.15,323.15,348.15,373.15 

2 

Hexane - Toluene - TMS 

298.15 

3 

Hexane - Xylene -TMS 

298.15 

4 

Octane - Banzene - TMS 

298.15 

5 

Octane - Toluene - TMS 

298.15 

6 

Octane - Xylene - TMS 

298.15 

7 

Heptane - Toluene - TMS 

298.15,323.15,348.15,373.15 

8 

Heptane - Xylene - TMS 

298.15,323.15 

9 

Cyclohexane - Benzene - TMS 

298.15,323.15,348.15,373.15 

10 

Heptane - Toleune - PC 

298.15 

11 

Heptane - Xylene - PC 

298.15 

12 

Heptane - Benzene - NFM 

1 

311.15 

13 

Heptane - T oluene - NFM 

311.15 

14 

Heptane - Xylene - NFM 

311.15 

15 

Flexane - Benzene - DMSO 

298.15,323.15 

16 

Heptane - Toluene - DMSO 

298.15 

17 

Hexane - Benzene - TREG 

293.15,313.15,333.15 

18 

Heptane - Toluene - NMP 

298.15 

19 

Cyclohexane - Benzene - DMF 

298.15 

20 

Heptane - Toluene - DEG 

298.15 


Table 5.1(b) Multicomponent systems used for comparison 


System 

Number 

Multicomponent systems 

Temperature (K) 

21 

Hexane - Benzene - Xylene - TMS 

298.15 

22 

Hexane - Octane - Benzene - TMS 

298.15 

23 

Octane - Toluene - Xylene - TMS 

298.15 

24 

Hexane - Octane - Benzene - Toluene - TMS 

298.15 

25 

Hexane - Heptane - Toluene - Xylene - PC 

298.15 
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Therefore the parameters obtained by this method can not be used for systems other 
than the one from which the parameters were regressed. Another problem with this 
method is that of initial guess. Because the parameter estimation is a minimization 
problem with several local minima’s, different initial guesses lead to the different set of 
parameters (Prausnitz et. al., 1980 [57]). The above two problems are very common in 
parameter estimation, as it can be observed from the results obtained for the systems 
Hexane-Benzene-Sulfolane and Hexane-Benzene-Dimethylsulfoxide. 

5.1.1 Hexane-Benzene-TMS 

The UNIQUAC binary interaction parameters at 298.15 K for the Hexane- 
Benzene-TMS have been regressed using the ternary LLE data reported by (Chen et. 
al., 2000 [13]). Many different initial guesses of parameters have been tried. Parameters 
obtained for two different initial guesses are reported in Table 5.2. This Table also 
shows the objective function values and rmsd values. The initial guesses of parameters 
used for the regression have also been reported in Table 5.2. The first initial guess has 
been obtained from Aspen Plus Simulator, while the second from the minimization of 
activity difference objective function of equation (4.3). It can be observed from Table 

5.2 that although two different initial guesses lead to two different sets of parameters, 
both sets represent the LLE data fairly well. The rmsd values and objective function 
values are also more or less same for the two initial guesses. Figure 5.1 and 5.2 
compare the experimental and predicted compositions on ternary diagrams. It can be 
seen from these figures that the predicted tie line compare well with experimental tie 
lines for both the initial guesses. 

5.1.2 Hexane-Benzene-DMSO 

A similar behavior is observed for this system also. Table5.3 gives the 
parameters regressed for this system, using two different initial guesses. The first initial 
guess is obtained through Aspen Plus, while the other is obtained using minimization of 
activity difference objective function of equation (4.3). The LLE data used for 
regression have been reported by ( Cassell et. al., 1989 [1 1] ) at a temperature of 298.15 
K. For this system also the two different initial guesses of parameters lead to two 
different converged set of parameters. The rmsd values for both sets are comparable 
with each other, but the parameters obtained are quiet different. 
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Hexane - Benzene - TMS at 298.15 K 
Predicted using initial guess I (Table 5.2) 



Predicted binodal dureve Hexane 

© 0 Predicted tie line 
0—0 Experimental tie line 

Fig: 5.1 Comparison of Predicted LLE with separate Estimation with initial guess I 
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Hexane - Benzene - TMS at 298.15 K 
Predicted using initial guess II (Table 5.2) 
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Table 5.2 UNIQUAC parameters for Hexane-Benzene-TMS at 298.15 K 


S.N. 

Binary pair 

Initial 

Regressed 

Objective 

No. of tie 




guess 

parameters 

Function 

lines used 

rmsd 


/ 

./ 

a,. 


-InZ. 

m 


1 

Hexane 

Benzene 

-116.10 

-369.0 





Hexane 

TMS 

561.40 

708.46 





Benzene 

Hexane 

157.95 

529.90 





Benzene 

TMS 

-21.90 

-27.31 

-159.47 

10 

0.6016 


TMS 

Hexane 

67.84 

133.40 





TMS 

Benzene 

59.08 

-41.20 




1 

Hexane 

Benzene 

27.64 

80.89 





Hexane 

TMS 

712.51 

646.52 





Benzene 

Hexane 

10.26 

-32.41 





Benzene 

TMS 

255.49 

1.052 

-143.87 

10 

0.8236 


'FMS 

Hexane 

87.52 

64.44 





'FMS 

Benzene 

-64.86 

113.45 





'Fable 5.3 UNIQUAC parameters for Hexane-Benzene-DMSO at 298.15 K 


S.N. 

Binary pair 

ij 

Initial 

guess 

Regressed 

parameters 

Objective 

Function 

No. of tie 
lines used 

rmsd 


i 

./ 



-Inl 



15 

Hexane 

Benzene 

-116.10 

-116.11 





Hexane 

DMSO 

956.60 

813.14 





Benzene 

Flexane 

157.20 

157.96 





Benzene 

DMSO 

84.24 

29.611 

-118.13 

7 

0.4954 


DMSO 

Hexane 

36.42 

186.98 





DMSO 

Benzene 

29.44 

-81.47 




15 

Hexane 

Benzene 

-155.46 

-31.47 





Hexane 

DMSO 

799.49 

4352.75 





Benzene 

Hexane 

84.614 

28.95 





Benzene 

DMSO 

-7.313 

88.46 

-104.48 

7 

0.3264 


DMSO 

Hexane 

55.456 

22.32 





DMSO 

Benzene 

-11.82 

-21.97 

1 




Furthermore, the parameters for common pair Hexane-Benzene between two 
systems are much different as summarized in Table 5.2 and 5.3. Therefore, parameters 
obtained using only one set of ternary data are found to be specific to that particular set 
of ternary system. Hence, these binary interaction parameters cannot be used directly to 
predict multi - component LLE of other systems involving these binaries. The 
parameters reported in DECHEMA (Sorensen and Arlt, 1979 [71]) also have the same 
problem. 
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To improve on above-mentioned problems, the parameter estimation procedure 
in this thesis has been modified so that the parameters between the common pairs in 
dilleient systems should be same. The teclmique used here is that the parameters for the 
systems involving common pairs are regressed simultaneously, as discussed in next 
section. 

5.2 Simultaneous Estimation of parameters involving common pairs 

In the technique used in this thesis, the objective function of equation (4.12) is 
calculated for each ternary system, while keeping the parameters between common 
pairs same in all systems. The sum of the objective functions of all the systems is then 
minimized. For example the interaction parameters between the Hexane -TMS pair are 
kept same in all the systems involving this pair. 

5.2. T Systems with TMS 

UNIQUAC parameters for Hexane-Benzene-TMS(l), Hexane-Toluene- 
TMS(2) and Hexane-Xylene-TMS(3) systems have been estimated simultaneously. 
The ternary LLE data reported by (Chen et. al., 2000 [13]), at 298.15 K temperature 
have been used for regression. The objective function of the regression has been the 
sum of individual objective functions for three systems. To represent each ternary 
system, 6 binary interaction parameters are needed. Therefore for three ternary systems 
1 8 binary parameters are required. But the Hexane-TMS binary pair is common to all 
three systems. Therefore for these three ternary systems only 14 parameters are 
required. In Table 5.4 the underlined components represent, the common binary pairs 
for which the interaction parameters were have to be same in all three systems. The 
initial guess of parameters have been obtained using UNIFAC-LL model of Aspen 
Plus. The regressed parameters along with the initial guesses are reported in Table 5.4. 
This Table also gives the objective function values after the regression. The predicted 
mole fractions are compared with experimental mole fraction, using the rmsd defined 
by equation (4.13) and shown in Figure 5.3. The predicted compositions, compares well 
with experimental compositions, as it can be observed from small rmsd values. Using 
the same approach UNIFAC group interaction parameters have also been obtained. The 
mutual interaction parameters for the groups CHi, ACH, and ACCH3 are taken directly 
from the UNIFAC LLE group interaction parameters given by ( Magnussen et. al., 
1981 [40] ) and are 
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Table 5.4 Simultaneous UNIQUAC parameters for Hexane - Benzene - TMS, Hexane 


- Toluene - TMS and Hexane - Xylene - TMS at 298. 1 5 K 


S.N. 

Binary pair 

ij 

Initial 

guess 

Regressed 

parameters 

Objective 

Function 

No. of tie 
lines used 

rmsd 


i 

J 

ay 

a 

V 

-Inl 

m 


1 

Hexane 

Benzene 

-116.10 

-155.47 





Hexane 

TMS 

561.40 

799.48 





Benzene 

Hexane 

157.95 

84.61 





Benzene 

TMS 

-21.90 

-7.31 

-160.06 

10 

0.5698 


TMS 

Hexane 

67.84 

55.46 





TMS 

Benzene 

59.08 

-11.82 




2 

Hexane 

Toluene 

167.27 

309.32 





Hexane 

TMS 

561.40 

799.48 





Toluene 

Hexane 

-137.20 

-207.88 

i 




Toluene 

TMS 

47.06 

136.04 

-183.23 

10 

0.3069 


TMS 

Hexane 

67.84 

55.46 





'FMS 

Toluene 

79.69 

-3.556 




3 

Hexane 

Xylene 

272.432 

303.15 





Hexane 

TMS 

561.40 

799.48 





Xylene 

Hexane 

-212.11 

-196.06 





Xylene 

TMS 

98.42 

181.20 

-157.80 

10 

0.7353 


TMS 

Hexane 

67.84 

55.46 





TMS 

Xylene 

65.01 

-5.14 





Total Objective function value = (-lnZ-|) + (-lni 2 ) + (“lnZ 3 ) = -502.1 


kept fixed at the values shown in italic and interaction parameters of these groups with 
the TMS solvent have been estimated an given in Table 5.5. Breakdown of components 
of the three ternary systems, into UNIT AC groups is given in Table 5.6. The underlined 
values of the parameters have been directly taken from (Magnussen et. al., 1981 [40] ). 
The predictions using regressed UNIQUAC parameters are better than the UNIQUAC 

as can be observed from rmsd values. 

The UNIT AC group parameters obtained using Hexane-Benzene-TMS, 
Hexane-Toluene-TMS and Hexane-Xylene-TMS can be used for the prediction of LLE 
for other systems also having paraffin, aromatic and TMS. The UNIFAC parameters 
obtained above has been directly used to predict the liquid-liquid equilibria for the 
systems Octane-Benzene-TMS, Octane-Toluene-TMS and Octane-Xylene-TMS 
systems. The predictions using these parameters is then compared with experimental 
data. Table 5.7 gives the comparison between experimental and predicted compositions 
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Hexane - Benzene TMS at 323.15 K 
Predicted using simultaneousiy regressed parameters (Table 5.17) 



4 




in terms of rmsd values and objective function values. The predictions of LLE using 
these UNIf AC parameters compare well with experimental data. In the similar way the 
parameters for three ternary systems, Octane-Benzene-TMS, Octane-Toluene-TMS 
and Octane-Xylene-TMS have been regressed. The LLE data have been taken from ( 
Chen et. ah, 2000 [13] ). Elere also the parameters for the underlined binary pair 
(Octane-TMS and TMS-Octane) common for all three systems remain same during the 
regression procedure. Furthermore the parameters for the italicized pairs Benzene- 
TMS, Toluene-TMS and Xylene-TMS are kept fixed at their earlier regressed values as 
obtained in fables 5.3. These parameters are not specific to the ternary system from 


fable 5.5 UNIFAC group interaction parameters for TMS with other main groups 


Binary Group Interaction Parameters 

j 

i 

CH 3 

CH, 

ACH 

ACCH 3 

TMS 

Clh 

0 

0 

-114.80 

-115. 70 

665.45 

CHt 

0 

0 

-114.80 

-115.70 

665.45 

ACM 

156.50 

156.50 

0 

167.00 

36.30 

ACCH 3 

104.40 

104.40 

-146.80 

0 

236.58 

TMS 

52.39 

52.39 

54.33 

9.93 

0 


Table 5.6 Comparison of Predictions using UNIFAC with experimental data for the 
systems Hexane-Benzene-TMS, Hexane-Toluene-TMS and Hexane-Xylene-TMS 


S. N. 

Ternary System 

UNIFAC 

Group assignment 

Objective 

Function 

rmsd 


Hexane 

2 CH 3 , 4 CHi 



1 

Benzene 

6 ACH 

-134.33 

0.9550 


TMS 

1 TMS 




Hexane 

2 CH 3 , 4 CH 2 



2 

'foluene 

5 ACH, 1 ACCH 3 

-161.18 

0.6022 


TMS 

1 TMS 




Hexane 

2 CH 3 , 4 CH 2 



3 

Xylene 

4 ACH, 2 ACCH 3 

-146.74 

0.9414 


TMS 

1 TMS 
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Table 5.7 Comparison of Predictions using UNIFAC with experimental data for the 
systems Octane-Benzene-TMS, Octane-Toluene-TMS and Octane-Xylene-TMS 


S. N. 

Ternary System 

Group assignment 

Objective 

Function 

rmsd 


Octane 

2 CH 3 , 6 CH 2 



1 

Benzene 

6ACH 

-121.25 

1.012 


TMS 

1 TMS 




Octane 

2 CH 3, 6 CH 2 



2 

Toluene 

5 ACH, 1 ACCH 3 

-136.93 

0.5208 


'TMS 

1 TMS 




Octane 

2 CH 3 , 6 CH 2 



3 

Xylene 

4 ACH, 2 ACCH 3 

-126.69 

0.7735 


TMS 

1 TMS 




where they have been regressed. Table 5.8 shows the regressed parameters at 298.15 K 
along with the initial guesses. The parameters obtained for systems (1,2 and 3 given in 
Table 5.4) having Hexane has been taken as initial guess. The predictions using 
regressed parameters compares well with experimental data as can be observed from 
reported rmsd. 

The validity of parameters regressed can be checked by comparison of 
prediction of multi-component LLE using these parameters with experimental data. 
Multi-component data for the quaternary systems Hexane-Benzene-Xylene-TMS and 
Hexane-Octane-Benzene-TMS have been used, because these systems involve same 
binary pairs for which the interaction parameters have been obtained. Experimental 
LLE data for these two systems have been reported by (Chen et. ah, 2000 [14]). Table 
5.9 compares the experimental compositions with predicted compositions using the 
IINIQIJAC parameters reported in Table 5.4. No parameters other than the reported in 
Table 5.4, are required for this quaternary system. The interaction parameters between 
Benzene and Xylene are kept zero, because this pair represents a fairly ideal solution ( 
Chen et. ah, 2000 [14] )■ Table 5.9 also gives the prediction of compositions using the 
UNIFAC method with regressed parameters obtained in Table 5.4. The ms'J values 
have been estimated for quaternary predictions and are reported in Table 5.9. For a 

quaternary system 
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rmsd = 100 


|l/2 


(5.1) 


w 4 2 

EEI 


k=\ /=i y=i 


8m 


Predictions of quaternary LLE using the regressed parameters have also been 
compared with experimental data for the system Hexane-Octane-Benzene-TMS. 
Experimental data have been reported by (Chen et. al., 2000 [14]). For prediction 
UNIQUAC parameters from Table 5.4 and Table 5.8 have been used. The interaction 
parameters between pair Hexane-Octane have been kept zero, as suggested by (Salem 
et. al., 1994 [70]). The predicted as well as experimental mole fractions are given in 
Table 5.10. These parameters have also been used to predict LLE for system Octane- 
Toluene-Xylene-TMS at 298.15 K for which the experimental LLE data have been 


Table 5.8 Simultaneous UNIQUAC parameters for Octane - Benzene - TMS, Octane - 


Toluene -- TMS, and Octane - Xylene - TMS at 298.15 K 


S.N. 

Binary pair 

Initial 

guess 

Regressed 

parameters 

Objective 

Function 

No. of tie 
lines used 

rmsd 


i 

j 



-Inl 

m 


1 

Octane 

Benzene 

-155.47 

-165.43 





Octane 

TMS 

799.48 

859.29 





Benzene 

Octane 

84.61 

111.67 





Benzene 

TMS 

-7.31 

-7.31 

-164.33 

10 

0.6485 


TMS 

Octane 

55.46 

34.23 





TMS 

Benzene 

-11.82 

-11.82 




2 

Octane 

Toluene 

309.32 

269.86 





Octane 

TMS 

799.48 

859.29 





Toluene 

Octane 

-207.88 

-190.21 





'Foluene 

TMS 

136.04 

136.04 

-158.95 

9 

0.3973 


TMS 

Octane 

55.46 

34.23 





TMS 

Toluene 

-3.56 

-3.55624 




.1 

Octane 

Xylene 

303.15 

320.43 





Octane 

TMS 

799.48 

859.29 





Xylene 

Octane 

-196.06 

-204.48 





Xylene 

TMS 

181.20 

181.20 

-154.54 

10 

0.5172 


TMS 

Octane 

55.456 

34.23 





TMS 

Xylene 

-5.14 

-5.14 





'Fotai Objective function value = (- In I, ) + (- In Zj ) + (“ In Z3 ) - 477.82 


AA 




reported (Chen et. al., 2000 [14]). The predicted mole fractions using the 
UNIQUAC and UNIFAC are given in Table 5.11 along with the experimental LLE 
data. The rmsd values are also given. 


Table 5.9 Comparison of experimental data for quaternary system Hexane(l)- 
Benzene(2)-Xylene(3)-TMS(4) with UNIQUAC and UNIFAC predictions. 


Tie 

line 


Raffinate phase 

Extract phase 




x[ 


xf 


x" 

1 

Exptl.* 

0.857 

0.076 

0.051 

0.017 

0.016 

0.052 

0.01 

0.921 

UNIQUAC 

0.874 

0.077 

0.047 

0.001 

0.017 

0.051 

0.015 

0.917 

UNIFAC 

0.869 

0.079 

0.049 

0.002 

0.020 

0.049 

0.012 

0.919 

2 

3 

I^xptl. 

0.709 

0.135 

0.143 

0.013 

0.021 

0.094 

0.043 

0.841 

UNIQUAC 

0.713 

0.139 

0.142 

0.006 

0.021 

0.091 

0.044 

0.844 

UNIFAC 

Fxpll. 

0.709 

0.138 

0.145 

0.008 

0.023 

0.091 

0.041 

0.845 

0.575 

0.266 

0.134 

0.025 

0.028 

0.199 

0.046 

0.727 

UNIQUAC 

0.576 

0.270 

0.135 

0.019 

0.029 

0.195 

0.045 

0.730 

UNIFAC 

0.577 

0.271 

0.133 

0.019 

0.031 

0.195 

0.047 

0.726 

4 

Exptl. 

0.450 

0.156 

0.356 

0.038 

0.026 

0.107 

0.135 

0.732 

UNIQUAC 

0.451 

0.154 

0.358 

0.037 

0.029 

0.109 

0.135 

0.726 

UNIFAC 

0.448 

0.154 

0.358 

0.040 

0.028 

0.109 

0.133 

0.729 

5 

ExptF 

0.316 

0.497 

0.098 

0.088 

0.083 

0.432 

0.062 

0.423 

“UNt^JAC' 

0.310 

0.493 

0.102 

0.095 

0.074 

0.432 

0.055 

0.439 

UNIFAC 

0.322 

0.495 

0.098 

0.085 

0.077 

0.434 

0.062 

0.427 

rmsd value using UN 

rmsd value using UNI 

rmsd value reported fc 

[QUAC = 0.588 

FAC = 0.4156 

)r NRTL = 0.69 (*Chen et. al., 2000 [14] ) 


Predictions using the regressed UNIQUAC and UNIFAC interaction 
parameters have been compared with experimental data for a five component system. 
LLE data for the quinary system Hexane-Octane-Toluene-Xylene-TMS have been 
reported by ( Chen et. al., 2000 [14] ). The predicted LLE data by UNIQUAC and 
UNIFAC are compared with experimental data in the Table 5.12. For prediction using 
UNIQUAC the parameters reported in Table5.4 and Table 5.8 have been used. The 
interaction parameters between psir Hexane-Octane and between pair Benzene- 
Toluene have been kept zero, as earlier. For prediction using UNIFAC the group 
parameters reported in Table 5.5 have been used. The vinsd vsXw&s are also given for 


comparison with reported rmsd values using NRTL specific parameters. Predictions of 
both UNIF AC and UNIQUAC are in well agreement with experimental LLE data. 

Using the similai’ approach the UNIQUAC parameters at 298.15 K for 
systems Heptane-Toluene-TMS and Heptane-Xylene-TMS have been estimated 
using the LLE data of Chen et al, 2000. The parameters for Toluene-TMS and Xylene- 
TMS pairs have been fixed at the values obtained earlier Table 5.3. Table 5.12 gives the 
regressed parameters along with initial guesses used. The rmsd values for the two 
systems are also reported in the Table 5.12. The parameters obtained for Hexane- 
Toluene-TMS and Hexane-Xylene-TMS have been used as initial guess. 


'Fable 5.10 Comparison of experimental data for quaternary system Hexane(l)- 
Octane(2)-Benzene(3)-TMS(4) with UNIQUAC and UNIF AC predictions. 


Tic 

line 


Raffinate phase 

Extract phase 


^2 

X3 

x' 

Xy 


xf 

x" 

1 

lixptl.* 

0.173 

0.655 

0.153 

0.019 

0.020 

0.007 

0.104 

0.869 

UNIQUAC 

0.189 

0.656 

0.153 

0.001 

0.004 

0.007 

0.104 

0.885 

UNIFAC 

0.186 


0.162 

0.002 




0.895 

2 


0.552 



0.021 

0.017 

0.003 

0.181 

mn 




ijgimi 

0.004 

0.018 

0.003 

0.179 






0.005 

0.019 

0.003 

0.167 

iggil 

3 

Exptl. 

0.288 


0.340 


0.011 

0.009 

0.260 

0.720 

UNIQUAC 

0.289 

0.357 

0.344 

Ijjgg 

0.013 



Iggg 

UNIFAC 

0.283 

0.351 

0.355 

Ijgig 

0.014 

ggyl 

0.244 



Exptl. 



llim 





gggg 

UNIQUAC 

IggM 




Ijgi^ 


iggi 


UNIFAC 

0.273 




|ig|H 

0.012 

Hggi 


5 

Exptl. 

0.154 

0.138 

0.618 



0.036 

gggi 

0.347 

UNIQUAC 

0.149 

0.130 

0.618 

0.102 

0.053 

0.039 

0.559 

0.349 

UNIFAC 

0.158 

0.141 

0.622 

0.079 

0.046 

0.030 

0.558 

0.366 

rmsd using UNIQl 

rmsd using UNIFA 

rmsd reported usir 

JAC = 0.7621 

.C = 1.05655 

ig NRTL = 0.99 (*Chen et. al., 2000) 


The LLE for the systems Heptane-Toluene-TMS and Heptane-Xylene-TMS have been 
pfedicted using the estimated UNIFAC parameters. The predictions have been 
compared with experimental data in Table 5.14. The rmsd values are comparable 
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with the nnsdvaluQS obtained using the UNIQUAC parameters regressed from the 
same set of data. These and previous results show the wide range of applicability of the 
regressed UNIT AC parameters. 

The same procedure is adopted for higher temperature data. The liquid-liquid 
equilibrium data given by Cassell et al, 1989 for Heptane-Toluene-TMS at 323.15 K 
and by Sorensen and Arlt, 1979 for Heptane-Xylene-TMS at 323.15 K have been used 
to estimate the UNIQUAC binary interaction parameters. The regressed parameters 
along with initial guesses are reported in Table 5.15. Parameters between underlined 
pairs are common to both the systems. Parameters at 298.15 K, reported in Table 5.13 
have been used as initial guess. The rmsd values and the objective functions are also 
reported in Table5.15. 


Table5.11 Comparison of experimental data for quaternary system Octane(l)- 
Toluene(2)-Xylene(3)-TMS (4) with UNIQUAC and UNIFAC predictions at 298.15 K. 


'I'ic 

line 


Raffinate phase 

Extract phase 

xi 

x[ 


x[ 

xf 

xf 

xf 


1 

Exptl.* 

0.649 

0.154 

0.184 

0.013 

0.007 

0.058 

0.041 

0.894 

UNIQUAC 

0.668 

0.155 

0.172 

0.005 

0.01 

0.06 

0.057 

0.873 






0.008 


mm 


2 

ExptL 

0.548 

0.244 

0.187 

0.022 


0.105 

0.057 

i^mil 

UNIQUAC 

0.559 

mnn 

0.183 






mmi 

0.82 

UNIFAC 

0.552 




0.009 

mm 

mm 




0.517 

0.182 



[ggHH 

mm 

0.088 

gmi 



mm 



gym 

mm 

mmi 

mm 

UNIFAC 

0.522 





mjm 

mmi 


4 

Exptl. 

^^1 




mm 

mm 


0.697 

UNIQUAC 

0.365 

0.422 



gym 

imH 

gmi 

0.689 

UNIFAC 

0.365 

0.415 



mm 


mmi 

0.689 

5 

Exptl 

0.251 

0.171 

1^^ 


mm 

imm 

mm 

immii 

UNIQUAC 

0.255 

0.17 

0.5 


mim 

immi 


0.67 

UNIFAC 

0.257 

0.166 



mm 


iimji 

0.668 

rmsd value using UN 
rmsd value using UNI 
rmsd value reported ft 

IQUAC = 0.81 

FAC = 0.62 

)rNRTL= 0.66 (* Chen et al., 2000 ) 
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Tables. 12 Comparison of experimental data for quinary system Hexane(l)-Octane(2)- 
Benzene(3)-Toluene(4)-TMS(5) with UNIQUAC and UNIFAC predictions. 


'fie line 

Raffinate phase 

Extract phase 

x; 

Xj 

Xj 

4 

4 

4 

xf 

xf 


xf 


Exptl.* 

0.696 

0.118 

0.083 

0.1 

0.004 

0.015 

0.002 

0.055 

0.042 

0.887 

1 

UNIQUAC 

0.695 

0.119 

0.084 

"o.ioi 

0.002 

0.016 

0.001 

0.054 

0.042 

0.887 


UNIFAC 

0.692 

0.119 

0.085 

0.101 

0.003 

0.018 

0.001 

0.053 

0.041 

0.886 


Exptl. 

L_ - 

0.232 

0.488 

0.17 

0.098 

0.012 

0.007 

0.009 

0.121 

0.043 

0.821 

2 

UNIQUAC 

0.232 

0.488 

0.175 

0.101 

0.004 

0.007 

0.007 

0.116 

0.039 

0.83 


UNIFAC 

0.23 

0.487 

0.179 

0.099 

0.006 

0.008 

0.007 

0.112 

0.042 

0.832 


I'Xpll. 

0.213 

0.243 

0.396 

0.117 

0.031 

0.019 

0.014 

0.31 

0.067 

0.59 

3 

UNIQUAC 

'0.208” 

0.237 

0.396 

0.125 

0.035 

0.02 

0.015 

0.308 

0.058 

0.6 


“iJNiixU'" 

0.212 

0.244 

0.397 

0.115 

0.031 

0.021 

0.013 

0.309 

0.069 

0.589 


l-xptl. 

0.154 

0.13 

0.133 

0.494 

0.088 

0.032 

0.02 

0.112 

0.334 

0.502 

4 

unIquac 

0.153 

0.13 

0.136 

0.496 

0.085 

0.029 

0.017 

0.107 

0.327 

0.519 


UNII'AC 

0.159 

0.136 

0.133 

0.49 

0.082 

0.03 

0.016 

0.113 

0.341 

0.501 


rnhsci value using UNIQUAC 


0.46 


nnsd value using UNIFAC = 0.37 

nm^ value reported for NRTL = 0.42 (*Chen et. al, 2000 [14] ) 


'Fable 5.13 Simultaneous UNIQUAC parameters for Heptane-Toluene-TMS and 


Ueptnae-Xylene-TMS at 298. 15K 


S.N. 

Binary pair 

ij 

Initial 

guess 

Regressed 

parameters 

Objective 

Function 

No. of tie 
lines used 

rmsd 


i 

./ 



-InZ 

m 



Heptane 

'Foluene 

309.32 

303.29 





1 leptane 

TMS 

799.48 

641.14 





Toluene 

Fleptane 

-207.88 

-205.85 




7 

'Foluene 

TMS 

136.04 

136.04 

-184.12 

10 

0.4070 


I'MS 

Hentane 

55.46 

48.28 





TMS 

'Foluene 

-3.56 

■ -3.56 





Heptane 

Xylene 

303.15 

281.55 





Heotane 

TMS 

799.48 

641.14 





Xylene 

Heptane 

-196.06 

-186.97 




8 

Xylene 

TMS 

181.20 

181.20 

-143.57 

10 

0.961 1 


TMS 

Fleotane 

55.46 

48.78 





TMS 

Xylene 

-5.14 

-5.14 





i'otal Objective function value s (-lnZ,) + (-lnZ, 2 ) = -327.69 
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Table 5.14 Comparison of Predictions using UNIFAC with experimental data for the 
systems Heptane-Toluene-TMS and Heptane-Xylene-TMS 


S. N. 

'Ternary System 

Group assignment 

Objective 

Function 

rmsd 


Heptane 

2 CHs, 5 CH 2 



7 

Toluene 

5 ACH, 1 ACCH 3 

-134.84 

0.8450 


TMS 

ITMS 



1 

Heptane 

2 CH 3 , 5 CH 2 



8 

Xylene 

4 ACH, 2 ACCH 3 

-125.53 

1.1516 


I’MS 

ITMS 




'Fable 5.15 Simultaneous UNIQUAC parameters for Heptane-Toluene-TMS and 


Ilcptnac-Xylcne-TMS at 323. 15K 



Binary pair 

Initial 

Regressed 

Objective 

No. of tie 


S.N. 

«.i 


guess 

parameters 

Function 

lines used 

rmsd 


i 

./ 



-Inl 

m 



1 Icptane 

Toluene 

303.29 

283.95 





llcptane 

TMS 

641.14 

600.12 





'Folucne 

Fleptane 

-205.85 

-210.69 




7 

'Foluene 

TMS 

136.04 

144.56 

-141.18 

10 

1.0160 


'FMS 

Heptane 

48.28 

40.23 





'i'MS 

Toluene 

-3.56 

-9.31 





Heptane 

Xylene 

281.55 

189.80 





Heptane 

TMS 

641.14 

600.12 





Xylene 

Heptane 

-186.97 

-161.93 




8 

Xylene 

TMS 

181.20 

165.98 

-77.87 

10 

1.5604 


'FMS 

Heptane 

48.78 

40.23 





'I'MS 

Xylene 

-5.14 

-2.80 




'Fotal Objective function value = (-lnL,) + (-lnI. 2 ) = ~219.06 


Liquid-liquid equilibrium data for system Heptane-Toluene-TMS at temperatures 
348.15 K and 373.15 K have been regressed. The LLE data reported by Sorensen and 
Arlt, 1 979 have been used for regression. For initial guess of parameters the parameters 
obtained at lower temperatures have been used. The parameters after regression are 
reported in Table 5.16. The rmsd values and objective function values are also reported 
in Table 5.16. 
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Heptane - Xylene - TMS at 323.15 K 
Predicted using simuitaneously regressed parameters (Table 5.15) 



Predicted binodal cureve 

©■ ■ •© Predicted tie line 
• Experimental tie line 


Fig: 5.4 Comparison of LLE Predictions for Heptane - Xyelne - TMS 
( Worst case, rmsd- 1.5604) 


0.0 


r A 



Table 5.16 UNIQUAC binary interaction parameters for Heptane-Toluene-TMS 


Temperature 

Binary pair 

ij 

Regressed 

parameters 

Objective 

Function 

No. of tie 
lines used 

rmsd 

i 

J 

a^j 

-Inl 

m 



Heptane 

Toluene 

303.29 




298.15 

Heptane 

TMS 

641.14 





Toluene 

Heptane 

-205.85 





Toluene 

TMS 

136.04 

-184.12 

10 

0.4070 


TMS 

Heptane 

48.28 





TMS 

Toluene 

-3.56 





Heptane 

Toluene 

283.95 



— 


Heptane 

TMS 

600.12 




323.15 

Toluene 

Heptane 

-210.69 





Toluene 

TMS 

144.56 

-141.18 

10 

1.0160 


TMS 

Heptane 

40.23 





'I'MS 

Toluene 

-9.31 





Heptane 

Toluene 

248.54 





Heptane 

TMS 

510.47 





Toluene 

Heptane 

-220.01 




348.15 

Toluene 

TMS 

164.52 

-147.85 

10 

0.8506 


TMS 

Heptane 

29.68 





TMS 

Toluene 

-33.73 





Heptane 

Toluene 

240.47 





Heptane 

TMS 

501.53 





Toluene 

Heptane 

-205.16 




373.15 

Toluene 

TMS 

98.96 

-131.44 

8 

0.7712 


TMS 

Heptane 

26.38 





TMS 

Toluene 

25.99 





Adopting a similar procedure UNIQUAC interaction parameters have been 
regressed for the system Hexane-Benzene-TMS at higher temperatures have been 
estimated. The experimental data have been reported by ( Sorensen and Arlt, 1979 
1 72]). Table 5.17 gives the regressed parameters for all the four temperatures 298, 
323.15, 348.15 and 373.15 K along with the values and objective function values. 

Using the similar approach the UNIQUAC interaction parameters for system 
Cyclohexanc-Benzene-TMS at different temperatures have been regressed. The LLE 
data of (Sorensen and Arlt, 1979 [72]) have been used. The regressed parameters axe 
given in Table 5.18. Parameters for pairs Benzene-TMS at different temperatures were 
kept fixed at their respective values obtained for system Hexane-Benzene-TMS. In 
I'able 5.18 the parameter values shown in italic were kept fixed at their respective 
values obtained for system Hexane-Benzene-TMS at different fehiperaWe^‘ Table 5.18 




also gives the objective function values at the regressed parameters along with the 
rwAJ values. For initial guess the parameters obtained for Hexane-Benzene-TMS 
system, given in Table 5.4 and Table 5.17 have been used. 


Table 5.17 UNIQIJAC binary interaction parameters for Hexane-Benzene-TMS 


Temperature 

Binary pair 

ij 

Regressed 

parameters 

Objective 

Function 

No. of tie 
lines used 

rmsd 

i 

J 


-InZ, 

m 



Hexane 

Benzene 

-155.47 





Hexane 

TMS 

799.48 




298. 15K 

Benzene 

Hexane 

84.61 





Benzene 

TMS 

-7.31 

-160.06 

10 

0.5698 


TMS 

Hexane 

55.46 





TMS 

Benzene 

-11.82 





I lexane 

Benzene 

-203.52 





1 lexane 

TMS 

715.93 





Benzene 

Hexane 

123.67 




323.15 K 

Benzene 

TMS 

-24.63 

-161.49 

10 

0.3312 


I’MS 

Hexane 

60.85 





TMS 

Benzene 

2.24 





Hexane 

Benzene 

-243.13 





Hexane 

TMS 

635.24 





Benzene 

Hexane 

193.70 




348.15 K 

Benzene 

TMS 

-24.51 

-134.95 

9 

1.3614 


TMS 

Hexane 

63.13 





TMS 

Benzene 

2.73 





Flexane 

Benzene 

-295.09 





Hexane 

TMS 

542.49 





Benzene 

Hexane 

229.44 




373.15 K 

Benzene 

TMS 

-24.24 

-150.62 

11 

1.2024 


TMS 

Hexane 

67.83 





TMS 

Benzene 

2.89 





5.2,2 Systems with PC 

As shown for the systems involving TMS as solvent that the parameters 
regressed simultaneously for three-four or more systems can be used for other systems 
also. Using the same approach as used for TMS, the binary parameters for systems 
involving PC have been regressed. These values of parameters obtained then used for 
prediction of multi-component LLE to compare with the experimental data. The 
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parameters for pairs that are common to systems with IMS as solvent have been kept 
same. Table 5. 1 9 gives the UNIQUAC interaction parameters regressed simultaneously 
Tabic 5.18 UNIQUAC paranieters for Cyclohexane-Benzeiie-TMS 


Temperature 

Binary pair 

ij 

Regressed 

parameters 

Objective 

Function 

No. of 

tie 

lines 

rmsd 

i 

j 

ay 

-]nL 

m 



Cyclohexane 

Benzene 

-328.34 





Cyclohexane 

TMS 

617.54 





Benzene 

Cyclohexane 

191.26 




298.15 K 

Benzene 

TMS 

-7.32 

-166.01 

10 

0.689 


TMS 

Cyclohexane 

70.22 





TMS 

Benzene 

-11.82 





Cyclohexane 

Benzene 

-263.83 





Cyclohexane 

TMS 

593.09 





Benzene 

Cyclohexane 

213.3157 




323.15 K 

Benzene 

TMS 

-24.63 

-168.40 

10 

0.455 


TMS 

Cyclohexane 

57.83 





TMS 

Benzene 

2.25 





Cyclohexane 

Benzene 

-296.61 





Cyclohexane 

TMS 

515.26 





Benzene 

Cyclohexane 

245.33 


j 


348.15 K 

Benzene 

TMS 

-24.51 

-153.34 

9 

0.392 


TMS 

Cyclohexane 

54.90 





TMS 

Benzene 

2.73 





Cyclohexane 

Benzene 

-336.93 





Cyclohexane 

TMS 

449.69 





Benzene 

Cyclohexane 

257.02 




373.15 K 

Benzene 

TMS 

-24.69 

-121.81 

8 

0.923 


TMS 

Cyclohexane 

52.15 





TMS 

Benzene 

3.22 





for Heptane-Toluene-PC and Heptane-Xylene-PC at 298.15 K. The experimental 
LLE data of Salem et al, 1994 and Salem, 1993 have been used for regression. The 
interaction parameters between the pair Heptane-Toluene and between pair Heptane- 
Xylene have been kept at the values given in Table 5.13. The values that were kept 
constant are shown in italic in Table 5.19. Interaction parameters for the pairs shown 
underlined are common to both the systems. For initial guess the parameters obtained 
for systems Heptane-Toulene-TMS and Heptane-Xylene-TMS respectively have been 
used. 
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I able 5.19 Simultaneous UNIQUAC parameters for Heptane-Toluene-PC and 
Heptane-Xylene-PC at 298. 15K 


S.N. 

Binary pair 

Initial 

Regressed 

Objective 

No. of tie 




guess 

parameters 

Function 

lines used 

rmsd 


i 

./ 



-InZ 

m 



Heptane 

Toluene 

303.29 

303.29 





Heptane 

PC 

836.30 

656.58 





Toluene 

Heptane 

-205.85 

-205.85 




10 

Toluene 

PC 

271.79 

237.09 

-120.84 

6 

0.1818 


P_C 

Heptane 

12.64 

7.43 


i 

1 


PC 

'Foluene 

-57.29 

-68.52 


! 

1 



/ leptane 

Xylene 

281.5 

281.50 



i 


1 Icplane 

PC 

836.3 

656.59 





Xylene 

Heptane 

-186.97 

-186.97 



I 

1 1 

Xylene 

PC 

633.54 

681.35 

-90.68 

6 

0.5888 


PC 

Heptane 

12.642 

7.43 





PC 

Xylene 

-128.82 

-160.79 





Total Objective function value = (-lni,) + (-lnZ,) = -211.52 


Table 5.20 UNIF AC group interaction parameters for PC with other main groups 


Binary Group Interaction Parameters a^ 


.j 

CH 3 

CH 2 

ACH 

ACCH 3 

PC 

CH3 

0 

0 

-114.8 

-115.7 

684.82 

CH 2 

0 

0 

-J14.8 

-115.7 

684.82 

AC II 

156.5 

156.5 

0 

167.0 

170.83 

ACCH 3 

104.4 

104.4 

! -146.8 

0 

525.39 

PC 

-8.42 

-8.42 

-36.91 

-136.62 

0 


rwW for system Heptane-Toluene-PC using above parameters = 0.4675 
rmsd for system Heptane-Xylene-PC using above parameters = 0.364806 


lire same systems as above have been used to estimate the UNIFAC group 
interaction parameters. Mutual interaction parameters between groups CH 3 , ACH, and 
ACCH 3 were directly taken from UNIFAC Table of (Magnussen et. al., 1981 [40]). 
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Heptane-Toluene-PC at 298.15 K 
Predicted using UlMiFAC (Table 5.20) 



Predicted Binodai curve 

© Predicted tie line 
• Experimental tie Hne 


^ 

Heptane 


Fig. 5. 5 LLE Predictions using regressed UNiFAC parameters 




The regressed group interaetion parameters are reported in the Table 5.20 along with 
the rmsd values. The underlined parameters were kept fixed at those values. 


I able 5.2 1 C ompaiison ol experimental LLE data with predicted for five component 
system Hexane ( 1 ) - Heptane (2) - Toluene (3) - Xylene (4) - PC (5) 


Tie 

- _ - Raftinate Phase mole Fractions 1 

Line 


y.1 


/ 

j 

/ 



X\ 

Xo 

X3 

X4 

X3 


Experimental 

0.482 

0.4472 

0.0284 

0.0405 

0.0019 

1 

Our Prediction 

0.4837 

0.446 

0.02817 

0.04093 

0.0012 


Reported Prediction 

0.5001 

0.4276 

0.0334 

0.0386 

0.0167 


Experimental 

0.456 

0.4152 

0.054 

0.0742 

0.0006 

2 

Our Prediction 

0.4567 

0.4132 

0.05316 

0.07517 

0.0018 


Reporlcd Prediction 

0.4660 

0.3985 

0.0623 

0.0726 

0.0006 


IvxperimenUil 

0.399 

0.366 

0.098 

0.136 

0.001 

3 

( )ur Prediction 

0.3977 

0.3607 

0.09945 

0.13851 

0.0036 


Reported Prediction 

0.4092 

0.3500 

0.1096 

0.1364 

0.0042 


I Experimental 

0.383 

0.3442 

0.1153 

0.1572 

0.0003 

4 

Our Prediction 

0.3815 

0.34 

0.11514 

0.15896 

0.0044 


Reported Prediction 

0.3853 

0.3297 

0.1292 

0.1572 

0.0045 


Experimental 

0.358 

0.324 

0.132 

0.185 

0.001 

5 

Our Prediction 

0.3543 

0.3185 

0.13267 

0.18878 

0.0057 


Reported Prediction 

0.3639 

0.3114 

0.1446 

0.1739 

0.0058 


Experimental 

0.3401 

0.3072 

0.1471 

0.2032 

0.0024 

6 

Our Prediction 

0.3376 

0.2997 

0.15215 

0.2036 

0.007 


Reported Prediction 

0.3408 

0.2917 

0.1650 

0.2032 

0.008 

Extract Phase mole fractions 



II 

11 

II 


II 



X, 

X2 

X3 

X, 

X5 


Experimental 

0.0189 

0.0081 

0.0099 

0.0115 

0.9516 

1 

Our Prediction 

0.0212 

0.0131 

0.0103 

0.0113 

0.9441 


Reported Prediction 

0.0107 

0.0121 

0.0053 

0.0123 

0.9719 


Experimental 

0.0189 

0.0081 

0.0188 

0.022 

0.9322 

2 

Our Prediction 

0.0213 

0.0129 

0.0199 

0.0214 

0.9245 


Reported Prediction 

0.0107 

0.0119 

0.0104 

0.0235 

0.914 


lExperimcntal 

0.0212 

0.0091 

0.0408 

0.0449 

0.884 

3 

Our Prediction 

0.0209 

0.0129 

0.0391 

0.042 

0.8851 


Reported Prediction 

0.0107 

0.0117 

0.0199 

0.0437 

0.8753 


1 Experimental 

0.0212 

0.0101 

0.0463 

0.0516 

0.8708 

4 

Our Prediction 

0.021 

0.0128 

0.0461 

0.0493 

0.8706 


Reported Prediction 

0.0107 

0.0116 

0.0244 

0.0528 

0.8643 


I Experimental 

0.0212 

0.0111 

0.0562 

0.0659 

0.8456 

5 

Our Prediction 

0.0208 

0.0129 

0.0546 

0.0605 

0.8512 


Reported Prediction 

0.0107 

0.0115 

0.0287 

0.0615 

0.8523 


i Experimental 

0.027 

0.0131 

0.0714 

0.0712 

0.8173 

6 

Our Prediction 

0.0208 

0.0128 

0.0639 

0.0671 

0.8353 


Reported Prediction 

0.0107 

0.0114 

0.0337 

0.0712 

0.8412 
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I ht itgrt/Sstd UNIFAC group interaction parameters have been used to predict 
a cjuinaiy system Hexane — Heptane — Toluene — Xylene — PC for which the 
experimental LLli data at 298.15 K have been reported by ( Salem et. ah, 1995 [69]). 

I he predicted mole fractions for this system using UNIFAC parameters reported in 
Table 5.20, are given in Table 5.21 along with experimental LLE data. Table 5.21 also 
gives the predicted LLE for this system by ( Salem et. al., 1995 [69] ) using specific 
UNIQUAC parameters. The predicted mole fractions are in good agreement with 
experimental data. The rmsd value for the predicted compositions reported by (Salem 
et. al., 1995 [69] ) is 0.876 as compared to 0.41 for this prediction. 

5.2.3 Sy.stenis with N-Formyl-Morpholine 

liinary interaction parameters for systems involving N-Formyl-Morpholine 
(NFM) as solvent have been estimated. Ternary experimental LLE data at 311.15 K for 
the systems I Icptane - Benzene - NFM, Heptane - Toluene - NFM and Heptane - 
Xylene - NMl*' have been reported by (Cinotti et. al., 1999 [16]). UNIQUAC binary 
interaction parameters for these three systems have been regressed simultaneously by 
adopting the similar approaches as used for TMS and PC systems. The regressed 
UNIQUAC parameters for these systems are reported in Table 5.23. Objective function 
values after the regression are also reported in this Table. The rmsd values reported in 
'I’able 5.23 for three systems, gives an indication of goodness of fit to LLE data. The 
initial guesses of parameters used are also reported in Table 5.23. 

The same ternary systems as used for UNIQUAC parameter estimation have 
been used to estimate the UNIFAC group interaction parameters. Group interaction 
parameters for NFM with three main groups CH3, ACH, and ACCH3 have been 
regressed. The mutual interaction parameters between these three groups have been 
kept fixed at there respective values given by (Magnussen et. al., 1981 [40] ). The 
regressed parameters along with the rmsdv&luss are reported in Table 5.22. The 
objective function values after regression for three systems have also been reported in 
'Fable 5.22. 1’he binary pairs common for the three ternary systems used are underlined 
in Table 5.23. The rmsd values obtained for UNIFAC are comparable with that of 
UNIQUAC. The total objective function value for the three ternary systems for 
UNIFAC’ is 1 93.06 as compared to -204.65 for UNIQUAC. 
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1 a.blC/ 5.22 UNIF AC gioup interaction parameters for NFM with other main groups 



Binary Group Interaction Parameters a,^ 

.i 

i 

CH 3 

CH 2 

ACH 

ACCH 3 

NFM 

CHa 

0 

0 

-114.8 

-115.7 

379.13 

CH 2 

0 

0 

-J14.8 

-115.7 

379.13 

ACH 

156.5 

156.5 

0 

167.0 

-2.05 

ACCH 3 

T04A 

104.4 

-146.8 

0 

113.61 

NFM 

38.18 

38.18 

78.67 

35.17 

0 

nusci and objective function for system F 

rmsd and objective function for system F 

nnsd and objective I'unction for system F 

leptane-Benzene-NFM = 0.6367, -60.33 

(eptane-Toluene-NFM = 0.7258, -56.149 

[eptane-Xylene-NFM = 0.2854, -76.5832 


Table 5.23 Simultaneous UNIQUAC parameters for Heptene-Benzene-NFM, Heptane- 


Toluene-Nl'M and Hcptane-Xylene-NFM at 31 1.5 


S.N. 

Binary pair 

ij 

Initial 

guess 

Regressed 

parameters 

Objective 

Function 

No. of tie 
lines used 

rmsd 


i 

./ 

% 


-In I 

m 


I2 

1 leptane 

Benzene 

303.29 

111.59 





1 Icntane 

NFM 

656.59 

363.60 





Benzene 

Heptane 

-205.85 

-153.63 





Benzene 

NFM 

237.08 

-17.96 

-61.86 

4 

0.5399 


NFM 

I leptane 

7.43 

46.55 





NFM 

Benzene 

-68.52 

5.10 




1.3 

1 leptane 

Toluene 

303.29 

293.16 





1 leptane 

NFM 

656.59 

363.59 





Toluene 

Heptane 

-205.85 

-208.34 





Toluene 

NFM 

237,08 

313.37 

-63.06 

4 

0.5426 


NFM 

I leptane 

7.43 

46.55 





Nl'M 

Toluene 

-68.56 

-130.58 




“14 

I leptane 

Xylene 

303.29 

233.84 





1 leptane 

Nl'M 

656.59 

363.59 





Xylene 

1 leptane 

-205.85 

-173.54 





Xylene 

NFM 

237.08 

112.94 

-79.70 

4 

0.2362 


Nl'M. 

1 leptane 

7,43 

46.55 





NFM 

Xylene 

-68.52 

8.36 





Total Obje 

ctive function value s 

(-lnX,) + (-lnT 2 ) + (~l’^ A) 

= -204.65 
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5.2.4 Systems with Dimethyl-Sulfoxide (DMSO) 

1 ernary experimental LLE for system Hexane - Benzene - DMSO at three 
different temperatures have been reported by (Cassell et. al., 1989 [11]). The LLE data 
at 298.15 K lor Hexane - Benzene - DMSO and Heptane - Toluene - DMSO have 
been used to estimate the UNIQUAC binary interaction parameters. Ternary LLE data 
toi system Heptane — Toluene — DMSO have been reported by (Sorensen and Arlt, 
1979 [72]). The estimated UNIQUAC parameters for Hexane - Benzene - DMSO at 
298.15 K and 323.15 K are given in Table 5.24 along with the initial guesses of 
parameters used. Ihe objective function and rmsd y&Iuqs are also reported in Table 
5.24. 'fhe rmsd values using the reported parameters by (Cassell et. al., 1989 [11]) are 
also given in Table 5.24. 


1 able 5.24 Ivstiniatcd UNIQUAC interaction parameters for system Hexane — Benzene 
- DMSO 



Binary pair 

Initial 

Regressed 

Objective 

No. of 


Temp. 



guess 

parameters 

Function 

tie lines 

rmsd 


/ 

./ 



-InZ, 

m 



Hexane 

Benzene 

-155.45 

-155.47 





Hexane 

DMSO 

544.64 

648.80 




298.15 

Benzene 

Hexane 

84.61 

84.61 





Benzene 

DMSO 

136.23 

129.69 

-112.05 

1 

0.6084 


DMSO 

Hexane 

58.75 

64.62 





DMSO 

Benzene 

-161.71 

-146.25 





Hexane 

Benzene 

-203.52 

-203.52 





I lexane 

DMSO 

542.42 

544.63 




323.15 

Benzene 

Hexane 

123.52 

123.55 





Benzene 

DMSO 

133.52 

136.23 

-114.19 

6 

0.4772 


DMSO 

1 lexane 

61.08 

58.75 





DMSO 

Benzene 

-151.57 

"•■■JT"' 

-161.71 





rm.v^/ using the UNIQUAC reported parameters at 298.15 - 1.578 


rmsd using the NRTL reported* parameters at 298.15 - 1.624 

rmsd using the UNIQUAC reported* parameters at 323.15 = 0.9472 
rmsd using the NRl'L reported* parameters at 323. 1 5 = 1 .0234 

* Reported by ( Cas.sell et. al., 1989 [1 1] ) 


UNIQUAC parameters for system Heptane — Toluene — DMSO at 298.15 K 
have been estimated and reported in Table 5.25. The rmsd vahxQS and objective 
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function values are also reported in Table 5.25 for this system. The interaction 
parameters shown in italic were kept fixed at the values, which have been regressed 
earlier. Regressed parameters tor Heptane — Toleuen- TMS system, were taken as 
initial guesses. 


Table 5.25 Estimated UNIQUAC parameters for Heptane - Toluene - DMSO at 298.15 


K. 


Temp. 

Binary pair 

ij 

Initial 

guess 

Regressed 

parameters 

Objective 

Function 

No. of 
tie lines 

rmsd 


i 

./ 


a,., 

-In I 

m 



Ilcpknw 

Toluene 

303.29 

303.29 





1 leplane 

DMSO 

641.14 

540.21 




298.15 

Toluene 

Ileplune 

-205.85 

-205.85 





Toluene 

DMSO 

136.04 

554.33 

-97.62 

8 

0.6075 


DMSO 

1 leptane 

48.78 

95.80 






DMSO 

'foluene 

-3.55 

-147.53 





Table 5.26 Simultaneous UNIFAC group interaction parameters for Hexane-Benzene- 
DMSO and 1 leptane-'roluene-DMSO at 298. 15K 


Binary (iroup Interaction Parameters a,, 


.i 

i 

CH3 

CH 2 

ACH 

ACCH 3 

DMSO 

CII 3 

0 

0 

-114.8 

-115.7 

722.11 

cih 

0 ^ 

0 

-114.8 

-115.7 

722.11 

At'll' 

] 56.5 

156.5 

0 

167.0 

153.99 

ACCIT,'" 

104.4 

104.4 

-146.8 

0 

188.46 


” 45.93 

45.93 

-44.48 

72.15 

T-'vx r\ i 

0 

00 nn Ar\AC 


nusd and objective function for system Hexane-Benzene-DMSO 0.61 83, -99.4046 
nnsd and objective function for system Heptane-Toluene-DMSO = 0.9452, -88.1820 


UNIFAC group binary interaction parameters for the groups CH 3 , ACH, 
and AC('l 1 i with DMSO have been regressed using the ternary LLE data for Hexane - 
Benzene - OMSO and Heptane - Toluene - DMSO. The regressed parameters are 
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given in I able 5.26. Mutual interaction parameters between the groups CH 3 , ACH, and 
ACCHj were kept fix at their respective values as given by (Magnussen et. al., 1981 
[401 )• fhc objective function and values are also reported in Table 5.26. The 
prediction using UNIFAC are comparable with UNIQUAC. 


5.2.5 Systems with Triethylene-Glycol (TREG) 

Liquid ~ Liquid equilibrium data for Hexane — Benzene — TREG have been 
reported by ( Sorensen and Arlt, 1979 [72] ). LLE data at 293.15 K, 313.15 K and 
3o.>.15 K temperatures have been reported. These data have been used to estimate the 
binary interaction parameters for the UNIQUAC model. For parameters estimation 
same procedure has been used, as used for other systems. The results of parameter 
estimations arc given in table 5.27. For initial guesses of parameters the activity 
difference objective function of equation (4.3) have been used. The LLE data at 293.15 
K have been used to minimize the activity difference objective function. The 
parameters obtained then were used as initial guess for objective function of equation 
(4.12). For parameter estimation at 313.15 K and 333.15 K the parameters obtained at 
lower temperatures have been used as initial guesses as shown in Table 5.27. The 
rmsd values and objective functions after regression are also reported in Table 5.27. 

5.2.6 Systems with N-Mcthyl-Pyrrolidone 

Ternary Liquid - Liquid equilibrium data for system Heptane - Toluene - NMP 

have been reported by (Ferreira et. al., 1984 [24]). LLE data at temperature 298.15 K 
have been used to regress the binary interaction parameters for UNIQUAC. The 
regressed parameters arc reported in Table 5.28. For initial guess of parameters the 
activity difference objective function of equation (4.3) have been minimized. The 
parameters obtained using activity difference objective function were then used to 
minimize the objective function of equation (4.12). The initial guess of parameters is 
also reported in Table 5.28. The values are also reported in Table 5.28 along with 
the objective function value. 

UNII'AC group interaction parameters for the main groups CH 3 , ACH, and 
ACCI I 3 with NMP have been estimated using the Heptane - Toluene - NMP at 298.15 
K. I.LF, Data rept)rted by (Ferreira et. al., 1984 [24]) have been used for this. The 
UNIFAC group interaction parameters between NMP and other main groups are 
reported in I'able 5.28. The mutual interaction parameters between the groups CH 3 , 
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ACH, and ACCH3 were kept fix at their respective values as given by Magnussen et aL, 
1981. 

5.2.7 Systems with Dimethyl-Formamide (DMF) 

1 ernary IJJ^, data lor the system cyclohexane - Benzene - DMF at 303.15 K 
have been reported by Sorenson and Arlt, 1979. UNIQUAC binary interaction 
pai-anieters for system Cyclohexane (CY) - Benzene - DMF have been estimated using 
these data. The estimated UNIQUAC parameters for this system are reported in Table 
5.29. The ttvisd values and objective function values are also given in Table 5.29. Due 
to unavailability of sufficient data for this system, UNIFAC parameters have not been 
estimated. 


I'ablc 5.27 I'.slimated UNIQUAC interaction parameters for system Hexane — Benzene 
- TRliC 



Binary pair 

Initial 

Regressed 

Objective 

No. of 


Temp. 

ij 


guess 

parameters 

Function 

tie lines 

rmsd 


/ 

j 

ay 

Qy 

-In I 

m 



1 lexane 

Benzene 

706.78 

262.23 





Hexane 

'I'REG 

805.26 

443.59 




293.15 

Benzene 

I lexane 

-299.73 

-161.23 





Benzene 

TREG 

346.03 

336.61 

-77.70 

4 

0.2818 


'FRliG 

Hexane 

68.68 

81.65 





TREG 

Benzene 

-103.28 

-77.48 





Hexane 

Benzene 

262.23 

233.92 





Hexane 

TREG 

443.59 

432.64 




313.15 

Benzene 

Hexane 

-161.23 

-145.64 





Benzene 

TREG 

336.61 

256.63 

-77.25 

4 

0.2690 


TRI'G 

Hexane 

81.65 

76.16 





'I'REG 

Benzene 

-77.48 

-36.33 





licxane 

Benzene 

233.92 






Hexane 

TREG 

432.64 

365.82 




333.15 

Benzene 

Hexane 

-145.64 

-183.29 





Benzene 

I’REG 

256.63 

246.01 

-74.98 

4 

0.3147 


'I'REG 

Hexane 

76.16 






TREG 

Benzene 

-36.33 

-29.17 
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Table 5.28 IJNIQUAC binary interaction parameters for the system Heptane - Toluene 
■ NMPat 298.15 K 


Temp. 

— .. 

Binar 

U 

y pair 

Initial 

guess 

Regressed 

parameters 

Objective 

Function 

No. of 
tie lines 

rmsd 

i 

./ 




-Ini 

m 



Heptane 

Toluene 

303.29 

303.29 





Heptane 

NMP 

641.14 

540.21 




298.15 

Toluene 

Heptane 

-205.85 

-205.85 





Toluene 

NMP 

136.04 

554.33 

-97.62 

1 

8 

0.6075 


NMP 

Heptane 

48.78 

95.80 





NMP 

Toluene 

-3.56 

-147.53 





Table 5.29 IJNIQUAC binary interaction parameters for the system Cyclohexane - 
Ben/enc DM h al 298,15 K 




1'emp. 

Binai 7 pair 

>j 

Regressed 

parameters 

Objective 

Function 

No. of 
tie lines 

rmsd 


/ 

./ 


-Ini 

m 



Cyclohexane 

Benzene 

354.79 





Cyclohexane 

DMF 

315.89 




298.15 

Benzene 

Heptane 

-228.99 





Benzene 

DMI' 

105.91 

-53.59 

3 

1.170 


DMl’ 

Cyclohexane 

3.23 



1 


DMF 

i 

Benzene 

-35.01 



i 

j 


Table 5.30 I INIFAC' group binary interaction parameters for NMP 


Binary Group Interaction Parameters 

J 

j 

CH3 

CH2 

ACH 

ACCH3 

NMP 

(TIj 

0 

0 

-114.8 

-115.7 

224.19 

(II 2 

0 

0 

-114.8 

-115.7 

224.19 

AC'li 

756.5 

156.5 

0 

167.0 

-130.68 

Af'Cllj'^”' 

104.4 

104.4 

-146.8 

0 

29.05 

NM.P 

19.78 

19.78 

76.52 

144.64 

O' 



rmsd = ' 

1.1003 




objective function = -1 54.8828 
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5.2.8 System with DEG 


Ternary UM data for the system Heptane - Toluene - DEG have been reported 
by Sorensen and Arlt, 1979. fhese data have been used to estimate the UNIQUAC 
binary parameters. Tiie regrc.ssed parameters are given in Table 5.31 along vvith initial 
guesses used. Parameters for the pair Heptane - Toluene were kept fix at the values 
obtained foi othei systems and shown italic in Table 5.31. The rmsd and objective 
function values are also reported in Table 5.31. 


I able 5.3 1 UNIQUAC binary interaction parameters for the system Heptane — Toluene 
.. .Dl'Xlat 298.15 K 


Temp. 

Binary pair 
ij 

Initial 

guess 

Regressed 

parameters 

Objective 

Function 

No. of 
tie lines 

rmsd 


i 

./ 



-Ini 

m 



1 Icptanc 

Toluene 

30129 

30129 





1 Icptanc 

DlXi 

641.14 

948.82 




298.15 

'foluenc 

1 Icptanc 

-205.85 

-205.85 





Toluene 

1)1X1 

136.04 

439.13 

-159.70 

8 

0.1699 


DlXi 

lieptane 

48.79 

-0.27 





DFXi 

Toluene 

-3.55 

-36.85 





5.3 Effect of Temperature on binary interaction parameters 

I'he effect of temperature on the interaction parameters for UNIQUAC has also 
been considered. To incorporate the effect of temperature on the binary interaction 
parameters, the parameters for UNIQUAC has been considered temperature dependent, 
as it has been used for UNIFAC parameters. The variation of parameters with 
temperature is shown in I-igurcs 5.6 to 5.10. The variation of parameters for Hexane - 
Benzene - TMS with temperature is shown in Figure 5.6. The parameters are reported 
in'I'ablc5.17. 

Similarly, the variation of parameters for systems Heptane - Toluene - TMS, 

Cyclohexane Benzene - TMS and Hexane - Benzene - DMSO with temperature is 

given in Figure 5.7, Figure 5.8, and Figure 5.9 respectively. The parameters are 
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reported in lable 5.16, I able 5.18, and Table 5.24 respectively. The temperature 
dependency ol UNIQUAC parameters considered here is given by equation (5.2) 

a,=a|;+a’(r) (4.21) 

T’hercl'ore a linear variation of parameters with temperature has been assumed. 
'The regression procedure has been modified to consider this temperature dependence of 
parameters. LLE data for the same systems at different temperature have been used. For 
example for system (1), the data at temperatures 298.15 K, 323.15 K, 348.15 K, and 
373.15 K have been used to regress simultaneously the temperature dependent 
parameters, 
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Interaction Parameter 


Variation of UNIQUAC parameters with temperature 



Temperature (K) 


Fig. 5.6 Variation of Binary Interaction Parameters with Temperature 
Hexane - Benzene - TMS 
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Interaction Parametere 


Variation of UNIQUAC parameters with TemperatLire 



Temperature (K) 

Fig: 5.7 Variation of Binary Interaction Parameters with Temperature 
Heptane - Toluene - TMS 


67 




Variation of UNIQUAC Parameters with Temperahire 



Fig. 5.8 Variation of binary interaction parameters with temperature 
For Cyclohexane - Bezene - TMS 
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Interaction Paranwter 


Variation of UNIQUAC parameters with Temperature 



Temperature (K) 


Fig:5.9 Variation of Binary Interaction Parameters with Tamperature 
Hexane - Benzene - DMSO 
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Chapter 6 

Solvent Performance Evaluation Results 

6.1 Interchange Energies 

As discussed in cluipler 3 the interchange energy given by equation (3.6) and 
(3.7), gives an indication ol solvent performance. These energies are estimated using 
the UNII'AC group parameters. I he estimated UNIFAC group parameters have been 
used lor this purpose, lable 6.1 gives the interchange energies calculated for five 
solvents using the parameters reported in chapter 5. 


Table 6. 1 Interchange Energies of Solvents with CHa and ACH groups 


Solvent 

Interchange Energies (K) 


CHz 

ACH 

TMS 

358.92 

45.32 

DMSO 

495.96 

56.95 

NME 

122.10 

-27.01 

Vi 

338.29 

66.93 

NFM 

208.65 

38.31 


It Ciin be observed from the 'Fable 6.1 that the difference in interchange energies 
of TMS with Cl I,’ and ACH groups, which represents nonaromatic and aromatic 
components respectively, is very large. Hence TMS will have a large selectivity 
towards arontatics. Similar behavior is observed for DMSO. The differences in 
interchange energies for other solvents are also appreciable, although lower than that of 
TMS and DMSO. 

6.1 Eerforiitaace of Solvents for a Single Stage Extraction 

Peii'ormance of different solvents for aromatic extraction has been evaluated. 
'I’he UNIQIIAC activity ct)eiricient model parameters obtained in previous chapter 
have been u.scd to estimate the selectivity, solvent power and performance index. To 
evaluate the performance of solvents the feed has been assumed to be of one aromatic 
key component and one non-aromatic key component. Selectivity and solvent power of 
solvents has been estimated for a single equilibrium extraction stage. For a particular 
feed composition, temperature and solvent to feed ratio the phase equilibrium 
compositions has been predicted using the UNIQUAC model with estimated 
liaramelers. For phase equilibria prediction the equations (4.16), (4.17) and (4.18) has 
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been .selveJ sinnillaneously. For this instead of taking mid point (z„ 2 ,) of tie lines, the 

feed composition has been used in equation (4.18), because the predicted tie line should 

pass through tlie feed point. 

Using tlic above procedure the selectivity, solvent capacity and performance 
indexes have been esUinated. for a feed of 20 % aromatics and for a solvent to feed 
ratio of 1.0, the selectivity and solvent capacity have been estimated at 298.15 K and 
are given in I able 6.2. Ihe aromatic - nonaromatic pair considered here is either the 
Hexane - Benzene or Heptane - Toluene, based on the availability of LLE data 

'l alde 6.2 Pcrforinance of Solvents for Single Stage Extraction 


Solvent 

SlccUvily 

Capacity 

PI 

TMS 

34.57 

0.711 

24.58 

DMSO 

22.25 

0.774 

17.24 

■PRPXi 

21.67 

0.533 

11.56 

NPM 

'12.44 

0.541 

6.73 

PU 

J2.92 

0.392 

5.07 

NMP 

3.34 

1.078 

3.69514 

DlXi 

7.22 

0.133 

0.967 


The variation ol‘ selectivity, .solvent capacity and performance index with aromatic 
content of Iced, solvent to feed ratio, and temperature have also been estimated. Figure 
6. 1 gives the variation of selectivity with aromatic content of the feed for different 
solvents for a solvent to feed ratio of 1.0 and temperature of 298.15 K. If has been 
observed (hat the selectivity for all solvents decreases with increasing aromatic content 
of feed. Selectivity ofTMS has been found greatest among the solvents considered 
here. Pure NMP, because of its small two phase region, cannot be used for high 
aromatic content feed. Figure 6.2 gives the variation of solvent capacity with aromatic 
content ol' feed, 'fhe solvent capacity of DMF is found to be greatest as compared to 
other solvents. Solvent capacity of NMP is higher than the other solvents for aromatic 
cotilcnt range ol' feed. The variation of performance index, based on the selectivity and 
solvent capacity, with aromatic content of feed is shown in Figure 6.3. 
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I he elttct of solvent to feed ratio on the solvent performance can be observed 
by considering a feed with 50 mole % aromatics at 298.15 K temperature. The effect of 
solvent to feed ratio on the selectivity of different solvents is shown in Figure 6.4.For 
almost all solvents the selectivity increases with solvent to feed ratio, except for DEG. 

1 he selectivity ol DEO does not change appreciably with changing solvent to feed 
ratio. The effect of solvent to feed ratio on solvent capacity is shown in Figure 6.5. The 
solvent capacity of NMP is found to be highest. The capacity does not ehange 
appreciably for almost all solvents. Similarly the variation of performance index with 
solvent to feed ratio is shown in Figure 6.6. The performance of each solvent increases 
with solvent to feed ratio, as may be expected. 

E fleet of temperature on performance of solvents has also been estimated. 
Idgure 6.7 gives the elTcct of temperature on selectivity of different solvents towards 
aromatics lor a feed containing 50 mole % aromatics and a solvent to feed ratio of 1.0. 
Selectivity of all tlic solvents decreases with temperature. Figure 6.8 provides the 
variation of performance index with temperature. The performance index decreases 
with temperature, except for PC, which shows a different behaviour then other 
solvents. Therefore the propylene carbonate seems have better extraction properties at 
higher temperatures. 





Selectivity 




Solvent Capacity 


Temperature = 298.15 K 
Solevnt to Feed Ratio = 1 



0 20 40 60 80 100 


% Aromatics in Feed 

Fig. 6.2 Variation of Solvent Capacity with % Aromatics in Feed 
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Fig; 6.4 Variation of Selectivity with Solvent to Feed Ratio 
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Fig:6.5 Variation of Solvent Capacity with Solvent to Feed Ratio 
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Fig;6.6 Variation of Performnce Index with T emperature 
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Fig; 6.8 Variation of Performance Index with Temperature 
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Chapter 7 

Conclusions and Recommendations 


Conclusions 

Methodologies have been established to estimate the binary interaction 
parameters from ternary liquid - liquid equilibrium data. Inside Variance Estimation 
melluKl based maximum likelihood function has been used for parameter estimation. 
'I'he UNlf'AC’ and DNIQUAC activity coefficient models have been considered. Both 
this models arc capable of reproducing liquid - liquid equilibria fairly well. It has been 
ohscrvctl that the dilferent set of initial guess of parameters may lead to different set of 
parameters. Although the dilTcrcnt set of parameters thus obtained for a particular 
system may preilict the liquid liquid equilibria for that system fairly well, that cannot 
be used for other systems. 

It htis been found that the simultaneous estimation of parameters improves on 
the above jiroblem. In llie parameter estimation technique used in this thesis, the 
parameters fu' common pair in different systems are estimated simultaneously. In this 
way the paramctei-s obtained are found to be generally applicable and not specific to the 
ternary system from where that have been regressed. The generality of parameters thus 
obtained has been tested for multicomponent systems. It has been shown that the 
parameters obtained simultaneously from different ternary systems can be used for 
prediction of multicomponent liquid - liquid equilibria. 

The lINUdUAC’ parameters estimated are found to be temperature dependent. 
The temperature dependency of parameters for the systems studied, have been 
considered. It has been observed that the temperature dependency behaviour of 
parameters is monotonic with linear relationship tor most of the cases. But for the some 
pairs the tem[>eraturc dependence does not follow the linear relation. 

Binary interaction parameters for the UNIQUAC and UNIFAC model from 
ternary liciuid liciuid equilibrium data have been estimated. Both these models can 
represent 1 ,1 .!■ fairly well for aromatic extraction systems. UNIFAC is found to be more 
suitable for multicomponent Iddi. 


m 1 



The extraction abilities of solvents used industrially and the solvents suggested 
by Molecular Design ot Solvents have been estimated. The solvents have been 

categorixed on the basis of a performance index, which incorporates both selectivity 
and solvent power. 


Recommendations 

An attempt should be made to discern the physical meaning of interaction 
parameters. Because these parameters are directly related with interaction energies 
between components or groups, this relationship along with experimental information 
can be usetl to get the binary parameters, which are not specific to ternary system used 
for their estimation. 

More systems should be incorporated into the simultaneous estimation 
techniciue, rather then to use only two or three systems. This may lead to more 
appropriate set of parameters, which can be used for multicomponent systems 
confidently. 

To incorporate temperature dependence in the parameters more proper initial 
guess can be obtained by extrapolating the parameter values obtained at other 
temperatures. 

UNIl'AC' parameters for many groups are not yet available for many groups. To 
use the molecular design of solvents technique efficiently, these parameters must be 
estimated. 

Water if present in small amount with solvents, is knoAvn to enhance the 
selectivity and to broaden the two phase region. Therefore effect of water on the 
performance of solvents must be considered. For a better performance, water can be 
mixed with solvents like NMP, which have a good capacity but poor selectivity for 
aromatics. 
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Appendix A 

Activity Coefficient Models 


A.l UNIQUAC Activity Coefficient Model 

The Universal Quasi Chemical (UNIQUAC) model, was developed by (Abrams 
& Prausnitz, 1975), to express the excess Gibbs free energy of mixtures. This model uses 
two adjustable parameters per binaiy. The UNIQUAC equation contains two parts, a 
combinatorial part and a residual part. The combinatorial part takes into account the 
composition, size and shape of the constituent molecules and contains pure component 
properties only. The residual part takes into account the intermolecular forces and 
contains two adjustable parameters, 'fhc molar excess Gibbs energy for a mixture is given 
by the sum of 


}''' (combinatorial) 

lu' 


/»! 


CD 


0. 


+ -y OjxAia- 

X, 2t;^' o, 


(A.1) 


and 


where 


g''' (residual) 
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cD; is segment or volume fraction of component i = 
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0: is area fraction of component / ^ 
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is volume parameter of component i; 
q, is surface area parameter of component i; 
T:, is adjustable parameter, given by 
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U.J is the interaction energy for the interaction of molecules of component i with 

the molecules of component / ; 

and Uy = Ujj 

a-j is the interaction parameter between component i and component j ; 


and Cy ^ Qy 

z is the coordination number, which is usually taken as 10. 

From equations (All. I ) and (A11.2) the activity coefficient for any component i , 
is given by 

In = In y[ (combinatorial) + In yf (residual) (A.4) 

where 


•n r / ' = In ^ + w/, In ^ ^ X x . 

X; 2 O,. 
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(A.6) 


'fhe pure component parameters are, respectively, measures of molecular van 
der Waals volumes and molecular surface areas. They are calculated as the sum of the 

group volume and group area parameters, and 

.,=y;vi'X; ,,=i;v<'>a (a-s) 

k k 

where , always an interger, is the number of groups of type k in molecule i. Group 
parameters Rj, and <2* normally obtained from van der Waals group volumes and 
surface areas as shown. 
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A.2 UNIFAC Activity Coefficient Model 

This group contribution method combines the solution of functional groups 
concept with a model for activity coefficients based on an extension of the quasi chemical 
theory of liquid mixtures, hence known as UNIFAC (UNIQUAC Functional-group 
Activity Coefficients). This model contains two adjustable parameters per pair of 
functional groups. The fundamental idea of a solution-of-groups model is to utilize 
existing phase equilibrium data for predicting phase equilibria of systems for which no 
experimental data are available. Like the UNIQUAC method, in this method also the 
activity coefficient of a component in a mixture is represented as product of a 
combinatorial and a residual activity coefficient. Hence the activity coefficient is given 
by 


In = In y. (combinatorial) + In yf (residual) (A.4) 

The combinatorial term is same as used for UNIFAC and given by (A.5). Only 
pure component parameters are needed for this part. Parameters r,. and are calculated 
from group parameters, and Q^. . These group parameters are obtained from the van 
der Waals group volume and surface areas and..^^ as given by (Reid et. al., 1988) 


R. 


n 


15.17 


and = 


A, 


2.5x10' 


The residual part, which is different from UNIQUAC and involves the 
interactions between groups, is given by 


inyf (a.9) 

k 


r* is the group residual activity coefficient and T^'^ is the residual activity 

coefficient of group A: in a reference solution containing only molecules of type i. The 
is the number of groups of type k in the molecule i . 

The group activity coefficient F^ is given by 
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Above equation (A. 10) also hold for InF^'^ . In this equation is the area 
fraction of group m , and the sums are over all different groups. ©„, is given by 



n 


(A.11) 


Where X,„ is the mole fraction of group in the mixture. The group interaction parameter 
'P,,,,, is given by 




mn 


= exp 


^mn ^nn 

RT 



(A. 12) 


Whs^"® U„,n is a measure of the energy of interaction bet'ween group’s m and n. The 

group interaction parameters must be evaluated from experimental phase equilibrium 
data. 
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Appendix B 

Simplex method for function minimization 


Nelder J. A. and R. Mead (Nelder and Mead, 1964) have described the simplex 
method for minimization of a function of « variables. Which depends on the 
comparison of function values at the («+l) vertices of a general simplex, followed by 
the replacement of the vertex with the highest value by another point. 

Consider the minimization of a function of n variables, without constraints. 

Pfi fhe («+l) points in « dimensional space defining the current 

“simplex”. Let the function value at is given hyy. , and define 

= m^(y,.) where suffix h is for “high” 


and 


y, = max(y;) where suffix / is for “low” 

/ 

Let the point P is the centroid of the points with j , and [PjPj} is the distance 

between P^ and Pj . At each stage in the process P^ is replaced by a new point; three 

operations are used - reflection, contraction, and expansion. These operations are 
defined as follows: 

The reflection of point P,, is denoted by P*, and its coordinated are defined by the 
relation 

Reflection : P* =(\. + a)P -aP^ 

Where a is a positive constant known as reflection coefficient. Thus the P* is on the 
line joining the points P;, and A on the far side of P from P^, with [P*P] = a[P^P]. If 
the function value at the reflection point that is y*, lies between y,, and y,, then 
P;, is replaced by P* and we start again with a new simplex. 

If y* < y, i.e. if reflection has produced a new minimum, then we expand 
P* to P** by the relation 
Expansion : p** _ ^p* ^ 

The expansion coefficient y , which is greater than unity, is the ratio of the distance 
[P**P] to [P*P]. 
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If y** <yi, we replace hy P** and restart the whole process; but if y** >>>/ 
then we have a failed expansion, and we replace by P* before restarting. 

If on reflecting P to P* we find thaty-* > y* for all i ^ h, i.e., replacing 
P by P* leaves y* the maximum, then we define a new P^ to be either the old P/^ or P* , 
whichever has the lower y value, and form 
Contraction: P** = pPf^ + (1 - /S)P 

The contraction coefficient yPlies between 0 and 1 and is the ratio of the distance 
[P**P] to [PP]. We then accept P** for P^ and restart, unless y** >min{y^,y*) , i.e. 
the contracted point is worse than the better of P^ and P* . For such a failed contraction 
we replace all the ’s by + P,)/2 and restart the process. 

Flow chart 

Calculate initial andy,. 
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Appendix C 

MATLAB Computer Programs 

A.l Program for initial guess using activity difference approach 

% FILE: initial_guess .m 

^icicicicic±icici(icicici(ici(icic-ki(i(iki(icicici(ikrk‘:kicicicicicicic‘ki(icic'k-:k-k'k'k-kic-ki(ic‘k‘ki(-k 

% INCLUDE FILE: act ivity_obj active .m 

% INCLUDE FILE: uniquac.m 

% INCLUDE FILE: exp_tie.txt 

% INCLUDE FILE: out . txt 

% SPECIFY temperature in this FILE 

% SPECIFY number of tie lines in FILE "activity_obj ective" 

% SPECIFY Rk Qk in FILE "activity_objective" 

^icicicicicicicicicicicicicicicic’kicicicic-kieicic'kificicificicicicicic'kic’k-k'kic'kicic-kic'k-kicicieicrk 

global exp_tie; 
global temp; 
global no_component ; 
global m; 

temp = 298.15; % temperature 
% OPEN FILE "exp_tie.txt 
fid =f open (' exp_tie . txt ') ; 

% READ tie line data in a MATRIX of (6,m) 
exp_tie = f scanf (f id, ' %g %g',[6 inf ] ) ; 

% CLOSE the FILE "exp_tie.txt" 
f close (fid) ; 

exp_tie = exp_tie'; % TRANSPOSE the MATRIX to (m,6) 

% INITIAL guess of parameters 
% ANY anitial guess will work 

a_interO = [263.884 798.177 -165.976 -49.8531 54.0035 176.911] 

} 

% Number of function evaluation 
options (14) =5000 ; 

% CALL FILE " act ivity_obj ective" for objective function 
% CALL "leastsq" for least square minimization 
a_inter = leastsq (' activity_obj ective a_inter0 , options) 

% OPEN FILE "out. txt" 

fiid = fopen ( ' out . txt ' , ' a ' ) ; % PRINT temperature 
fprintf (f iid, ' temprature = %g \n',temp); 

% PRINT initial guess 

fprintf (fiid, ' initial guess = %g %g %g %g %g %g\n' , a_inter0) ; 
% PRINT parameters 

fprintf (fiid, ' parameters_act = %g %g %g %g %g %g\n ' , a_inter) ; 
% CLOSE the FILE "out. txt" 
fclose (fiid) ; 
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^ icicificicicicicicicicicicicicicicic-kicic'k'k'k'k'kit-kic-k'kicrkic'kicicicic-k-kicicicrkic 

% FILE: uniquac.m 

% ‘k-kic-kick'k-kic'k'kicic-kic'k'kicieicic'kicic'kic-k'k'k-Jc'kicie'h-k-k'kieic'k'k-kiicie-k'k 

function gammal = gamma (no_component ,x, 
temp, Rk, Qk, a_interaction) 

% CALCULATE zi(i,j) 
for i = 1 : no_component 

for j = 1 : no_component 

z_i (i, j ) =exp (- (a_interaction (i, j ) /temp) ) ; 

end 

end 

sum_l = 0 ; 
sum_2 = 0 ; 

% SUMMATION x(i)*R(i) 

% SUMMATION x(i)*Q(i) 
for i = 1 : no_component 
theta_zi(i) = 0; 
sum ( i ) = 0 ; 

sum_l=sum_l+x (i) *Rk (i) ; 
sum_2 =sum_2 +x ( i ) *Qk ( i ) ; 

end 

sum_x_l = 0 ; 

for i= 1 : no_component 

phi (i) = (x (i) *Rk (i) ) /sum_l ; 
theta (i) = (x (i) *Qk (i) ) / sum_2 ; 

1 (i) =5* (Rk{i) -Qk(i) ) +1-Rk(i) ; 
sum_x_l = sum_x_l + x ( i ) * 1 ( i ) ; 

end 

for i=l :no_component 

for j =1 : no_component 

theta_zi (i) =theta_zi (i) + theta (j ) *z_i (j , i) ; 
end 

end 

% SUMMATION 

for i=l : no_component 

for j = 1 : no_component 

sum(i) =sum(i) + (theta (j ) *z_i (i , j )) /theta_zi (j ) ; 

end 

end 

% ESTIMATE gamma vector 
for i=l : no_component 

ln_gamma_r (i) = Qk(i)*(l-log( thet a_z i ( i ) ) - sum ( i ) ) ; 
ln_gamma_c (i) =log (phi (i) /x (i) ) +5*Qk (i) * (log (theta (i) /phi (i)))+ 
1 ( i ) - (phi ( i ) /x ( i ) ) * sum_x_l ; 

gammal (i) =exp (ln_gamma_c (i) +ln_gamma_r (i) ) ; 

end 
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%********************************************************** 
% FILE unifac.m 

%********************************************************** 
function gammal = unifac (no_component , no_group, v, x, 
temp, Rk, Qk, a) 

% INITIATE 
suml = 0 ; 
sum2 = 0 ; 
sums = 0 ; 
sum4 = 0 ; 
sum6 = 0 ; 

% no_compoent = Number of components 
% no_group = Total number of groups 
for i=l : no_component 
for k=l : no_group 
sum77 (i , k) =0 ; 
sum88 (i , k) =0 ; 

end 

end 

for i=l :no_component 
r(i) = 0; 
q ( i ) = 0 ; 

sumll (i) =0 ; 
sum66 (i) =0 ; 
sums (i) =0 ; 
sumlO (i) =0 ; 
end 

for j =1 : no_group 
sums ( j ) =0 ; 
sum7 ( j ) = 0 ; 
sum8 ( j ) =0 ; 

end 

for i=l : no_component 
for j =1 : no_group 

r(i) = r (i) +v(i, j) *Rk(j) ; 
q(i) = q{i) +v (i , j ) *Qk ( j ) ; 
sum4 = sum4 +v ( i , j ) *x ( i ) ; 

end 

end 

for i=l : no_component 

suml = suml+r (i) *x (i) ; 

sum2 = sum2 +q ( i ) *x ( i ) ; 

l(i) = 5* (r(i) -q(i) ) - (r(i) -1) ; 

sums = sumS+x (i) *1 (i) ; 

end 

for i=l : no_component 
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ln_gamma_c (i) = 

log(r (i) /suml) +5*q(i) *log( (g(i) *suml) / (r (i) *sum2)) . . . 

+1 (i) - (r (i) *sum3) /suml; 

end 

for j =1 :no_group 

for i=l : no_component 

sum5 (j ) =sum5 (j ) +v(i, j ) *x(i) ; 

end 

end 

for j =1 : no_group 

X ( j ) =sum5 ( j ) / sum4 ; 
sum6 = sum6+Qk ( j ) *X ( j ) ; 

end 

for j = 1 : no_gr oup 

xi (j) = (Qk(j) *X(j) ) /sum6; 

end 

for k=l : no_group 
for j =1 : no_group 

zi ( j , k) =exp (-a ( j ,k) /temp) ; 
sum7 (k) =sum7 (k) +xi ( j ) *zi ( j , k) ; 

end 

end 

for k=l;no_group 
for j =1 : no_group 

sums (k) = sums (k) + (xi ( j ) *zi (k, j ) ) /sum7 ( j ) ; 

end 

end 

for k= 1 : no_gr oup 

ln_gamma_k (k) = Qk (k) * (1-log (sum7 (k) ) -sumS (k) ) ; 

end 

for i=l :no_component 
for j = 1 : no_group 

sumll (i) =sumll (i) +v (i , j ) ; 

end 

end 

for i=l : no_component 
for j =1 : no_group 

XX (i, j ) =v (i, j ) / sumll (i) ; 
sum66 (i) =sum66 (i) +Qk (j ) *XX (i, j ) ; 

end 

end 

for i=l : no_component 
for j =1 : no_group 

xixi ( i , j ) =Qk ( j ) *XX ( i , j ) / sum6 6 ( i ) ; 
end 

end 

for i=l :no_component 
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for k=l : no_group 
for j =1 : no_group 

sum77 (i, k) =sum77 (i, k) +xixi (i , j ) *zi ( j , k) ; 

end 

end 

end 

for i=l :no_component 
for k=l:no_group 
for j =1 : no_group 

sums 8 (i,k) =sum88 (i, k) + (xixi (i,j)*zi(k,j)) /sum77 (i , j ) ; 

end 

end 

end 

for i=l :no_component 
for k=l:no_group 

ln_gamma_i_k (i,k) =Qk(k) * (1-log (sum77 (i,k) ) -sum88 (i,k) ) ; 

end 

end 

for i=l :no_component 
for j =1 : no_group 

sum9 ( i ) =sum9 ( i ) +v ( i , j ) * ln_gamma_k ( j ) ; 
suml 0 ( i ) =suml 0 ( i ) +v ( i , j ) * ln_gamma_i_k ( i , j ) ; 

end 

ln_gamma_r ( i ) = sum9 ( i ) - suml 0 ( i ) ; 

gammal ( i ) = exp ( ln_gamma_c ( i ) +ln_gamma_r ( i ) ) ; 

end 
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%************************************************************ 
% FILE: nrtl.m 

%******************-k-k-k*************-k*-k***i<-k******-k-k**-k**-k**** 

function gammal = gamman (no_component ,x, temp, alpha 

, a_interaction) 

for i = 1 : no_component 

for j = 1 : no_component 

tau (i, j ) =a_interaction (i, j ) /temp; 

G(i, j) =exp (- alpha* tau(i,j) ) ; 

end 

end 

for i=l :no_component 
suml ( i ) = 0 ; 
sum2 ( i ) = 0 ; 
sum3 ( i ) = 0 ; 

end 

for i=l :no_component 
for j = 1 : no_component 

suml (i) =suml (i) +tau (j , i) *G(j , i) *x(j ) ; 
sum2 (i) =sum2 (i) +G ( j , i) *x ( j ) ; 

end 

end 

for i=l :no_component 
for j =1 : no_component 

sum3 (i) =sum3 (i) + ( (x ( j ) *G (i, j ) ) /sum2 (j) ) * (tau(i, j) - 
suml ( j ) /sum2 ( j ) ) ; 
end 

end 

for i=l : no_component 

gammal (i) = exp ( (suml (i) /sum2 (i) ) +sum3 (i) ) ; 

end 
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A.2 Program for estimation of parameters from ternary LLE data for a 
particular ternary system 

%*******‘*'***-*‘********************:^**^:^********.************* 

% FILE: initial. m 

%********************************************************** 

% INCLUDE FILE: obj active. m 

% INCLUDE FILE: estimation . m 

% INCLUDE FILE: lle_equat ion . m 

% INCLUDE FILE: uniquac.m % same as for initial guess 

% INCLUDE FILE: unifac.m % same as above 

% INCLUDE FILE: nrtl.m % same as above 

% INCLUDE FILE: exp_tiel.txt 

% INCLUDE FILE: out . txt 

% SPECIFY temperature in "estimation" FILE 
% SPECIFY number of tie lines in "objective 
% SPECIFY Rk and Qk in estimation FILE 

global pi; % DECLARE global to objective function 

global rms; % DECLARE global to rmsd 

global itr; % DECLARE global to number of ITERATIONS 

itr =0; % Start ITERATION number 

% SET number of function evaluation 
options (14) = 5000; 

% Initial Guess of Parameters 
% Arrange the parameters like 

% a_inter0 = [ (1,2) (1,3) (2,1) (2,3) (3,1) (3,2)]; 

a_inter0 = [-116.105 590.329 157.958 13.5982 78.543 62.2065]; 

% CALL the "objective" FILE for objective function 

% CALL "fmins" for minimization 

a_inter = fmins (' obj ective ’, a_inter0 , options) 

% OPEN FILE "out. txt" for output 
fiid = fopen ( ' out . txt ' , ' a ' ) ; 

% PRINT initial guess used 

fprintf (fiid, ' initial guess = %g %g %g %g %g %g \n' , a_inter0) ; 
% PRINT parameters estimated 

fprintf (fiid, 'parameters = %g %g %g %g %g %g %\n ' , a_inter) ; 

% PRINT objective function 

fprintf (fiid, 'L = %g \n',pl) ; % PRINT rmsd values 
fprintf (fiid, ' pi = %g \n',rms); 

% PRINT iteration numbers 

fprintf (fiid, ' iterations = %d',itr); 

% CLOSE FILE "out. txt" 
f close (fiid) ; 
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^•k'k-k'k-k-k’k'k'kic-kicic'k'k'kic'k-kieic'k'k'kic-kicic-kickicic'k'kickicic-k'k'k'kicickicicicicick'kic 

% FILE: objective 

%*******j«r********************ilf*******-******* ************* 

function L = obj ective (a_inter) % FUNCTION DECLARE 
global pi ; 
global rms; 
global itr; 

m=10; % Number of experimental tie lines 

% INCREASE Number of iteraton by 1 
itr = itr+1 

% CONVERT the parameters in MATRIX form 
% a_interaction = [0 (1,2) (1,3) 

% (2,1) 0 (2,3) 

% (3,1) (3,2) 0 ]; 

a_interaction = [0 a_inter(l) a_inter(2); 

a_inter(3) 0 a_inter(4); 
a_inter(5) a_inter(6) 0]; 

% OPEN FILE "exp_tie.txt" for Experimental data 
fid = fopen ( ' exp_tie . txt ' , ' r ' ) ; 

% READ tie line data in (6,m) MATRIX 
exptiel = f scanf (f id, ' %g %g' , [6 inf ] ) ; 

% CLOSE "exp_tiel.txt" 
f close (fid) ; 

% TAKE the transpose of MATRIX to make it (m, 6) 
exp_tie = exptie'; 

% INITIATE sum's 

suml=0 ; sum2=0 ; sum3=0 ; sum4=0 ; 

% INITIATE v's 
vl=0 ; v2=0 ; V3=0 ; v4=0 ; 

% CALL "estimation" for tie line estimation 
% STORE them in a MATRIX 

est_tie = estimation (m, exp_tie , a_interact ion) ; 

% CALCULATE errors 
error = exp_tie - est_tie; 

% CALCULATE sum of errors 
for i = l:m 

suml= suml+error (i , 1) ; 
sura2= sum2+error (i, 2) ; 
sum3= sum3+error (i, 4) ; 
sum4= sum4+error (i , 5) ; 

end 

% AVERAGE of errors 

suml = suml/m; sum2 = sum2/m; sum3 = sum3/m; 
sum4 = sum4/m; 

% ESTIMATE variances 
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for i = l:m 

vl = vl+ (error (i , 1) -suml) '^2 ; 
v2 = v2+ (error (i , 2) -sum2) ^2 ; 
v3 = v3+ (error (i , 4) -sum3) ^2 ; 
v4 = v4+ (error (i, 5) -sum4) '^2 ; 

end 

vl=vl/ (m-1) ;v2=v2/ (m-1) ;v3=v3/ (m-1) ;v4=v4/ (m-1) ; 

% INITIATE 
tern = 0 ; 

% CALCULATE summation 
for i=l:m 

tern = tern + (error (i, 1) ) ^2/vl+ (error (i, 2) ) ^2/v2+. . . 

( error (i,4))^2/v3+( error ( i , 5 ) ) ^2 /v4 ; 

end 

% CALCULATE the objective function 

L = 2*m*log (2*pi) + (m/2) *log (vl*v2*v3*v4) + (1/2) *tem 
pi = L; 

sum_square = 0 ; 

% CALCULATE the rmsd 
for i=l:m 
for j=l:6 

sum_square = sum_square+ (est_tie (i , j ) -exp_tie (i , j ) ) ^2 ; 

end 

end 

rms = 100* (sum_square/ (6*m) ) ^0 . 5 
% RETURN 

% FILE: estimation. m 

^★★*****^*****‘3«r**-:*:*************:i<f*******‘*:«f********* 

% FUNCTION DECLARE 

function estimate = estimation! (m, exp_tie, a_interaction) 

% Rk for three components 
Rk = [4.4998 3.19051 4.0356]; 

% Qk for three components 
Qk = [3.856 2.4 3.2 ] ; 

% SPECIFY temperature 
temp =298.15; 
for i=l:m 

% ESTIMATE mid point of tie lines 

z_l = (exp_tie (i , 1) +exp_tie (i , 4) ) /2 ; 
z_2 = (exp_tie (i , 2) +exp_tie (i , 5) ) /2 ; 

% initial guess of a tie line 
XO = exp_tie (i, : ) ; 

% CALL "lle_equation" for setting LLE equations 
% CALL "fosolve" for Solving equations 
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X = 

f solve ( ' lle_equation ' , XO , [] , [] , z_l , z_2 , Rk, Qk, a_interaction, tern 
p) ; 

% STORE the tie lines in a MATRIX 
for j=l:6 

est (i, j) = X(j) ; 

end 

estimate = est; 

end 


^ icicicicicicicicicicicicicicicici(icicic9cici(i(icic’k‘kicic"k‘k-kic'ki(‘kici('ki('k-k‘kicicic'kicici(‘k'kic'kic 

% FILE: lle_equation.m 

^ ieicici(icicicicicici(icicicicic-kicici(i(i(icicicic-kici(icicrki(i(ici(ici(ic'kikicicicicic-:k'k'k'k'kic‘k'k’k'k 

% FUNCTION DECLARE 
function equation = 

lle_eqution (X, z_l, z_2 , Rk, Qk, a_interaction, temp) 

% number of components 
no_component = 3 ; 

% two phase compositions 

x{l) = X(l) ; 

x(2) = X{2) ; 

x{3) = X(3) ; 

y (1) = X{4) ; 

y (2) = X(5) ; 

y(3) = X(6) ; 

% CALL "uniquac" for activity coefficints for PHASE I 
gamma_I = uniquac (no_Gomponent , x , 
temp, Rk, Qk, a_interact ion) ; 

% CALL "uniquac" for activity coefficints for PHASE II 
gamma_II = uniquac (no_component , y , 
temp, Rk, Qk, a_interact ion) ; 

% SET the LLE equations 
equation = [ 

(z_l-X(l) ) * (X(5) -X(2) ) - (X(4) -X(l) ) * (z_2-X(2) ) ; 

X (1) *gamma_I (1) - X (4) *gamma_II (1) ; 

X(2) * gamma_I(2) - X(5) * gamma_II (2) ; 

X(3) * gamma_I(3) - X(6) * gamma_I I ( 3 ) ; 

X(l) + X(2) + X(3) -1; 

X(4) + X(5) + X(6) -1] ; 

% RETUREN 
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A.3 Program for Simultaneous estimation of parameters from three 
ternary systems 

% FILE: initial_simult .m 

^ikicicici(icicici(icicicici<ic‘kicificici(icicicicici(icicicicicicici(icicic-k-k-kic-k‘k‘k-k9cici('k‘k-k‘kicicicic‘k 

% INCLUDE FILE: obj ective_simult . m 

% INCLUDE FILE: estimation! . m 

% INCLUDE FILE: estimationC . m 

% INCLUDE FILE: estimations . m 

% INCLUDE FILE: lle_equation . m % same as for above 

% INCLUDE FILE: uniquac.m % same as for above 

% INCLUDE FILE: unifac.m % same as for above 

% INCLUDE FILE: nrtl.m % same as for bove 

% INCLUDE FILE: exp_tiel.txt 

% INCLUDE FILE: exp_tie2.txt 

% INCLUDE FILE: exp_tie3.txt 

% INCLUDE FILE: out . txt 

% SPECIFY temperature in all "estimation" FILE 

% SPECIFY number of tie lines in "obj ective_simult 

% SPECIFY Rk and Qk in each estimation FILE 

global pi; % DECLARE global to objective function 1 

global p2 ; % DECLARE global to objective function 2 

global p3 ; % DECLARE global to objective function 3 

global rmsl; % DECLARE global to rmsdl 

global rms2 ; % DECLARE global to rmsd2 

global rms3 ; % DECLARE global to rmsdS 

global itr; % DECLARE global to number of ITERATIONS 

itr =0; % Start ITERATION number 

% SET number of function evaluation 

options (14) = 5000; 

% Initial Guess of Parameters,,,,, 

% Dont Repeat parameters which are common 
% Arrange the parameters like 

% a_inter0 = [ (1,2) (1,3) (2,1) (2,3) (3,1) (3,2) % 

For First System 

% (1,2) (2,1) (2,3) (3,2) % 

For Second System 

% (1,2) (2,1) (2,3) (3,2)]; % 

For Third System 

% NOTE that parameters (1,3) and (3,1) are common to all 
systems 

a__inter0 =[-116.105 590.329 157.958 13.5982 78.543 62.2065... 
167.265 -137.203 80.06 37.5965... 

272.432 -212.1113 120.6137 56.2235]; 

% CALL the "obj ective_simult " FILE for objective function 
% CALL "fmins" for minimization 
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a_inter = f mins (' obj ective_simult a_interO , options) 

% OPEN FILE "out.txt" for output 
fiid = fopen ( ' out . txt ' , ' a ' ) ; 

% PRINT initial guess used 

fprintf (fiid, 'initial guess = %g %g %g %g %g %g %g %g %g %g %g 
%g %g %g\n' , a_interO) ; 

% PRINT parameters estimated 

fprintf (fiid, 'parameters = %g %g %g %g %g %g %g %g %g %g %g %g 
%g %g\n' ,a_inter) ; 

% PRINT objective functions 

fprintf (fiid, ' L = %g %g %g\n' ,pl,p2,p3) ; 

% PRINT rmsd values 

fprintf (fiid, ' pi = %g %g %g\n' ,rmsl,rms2,rms3) ; 

% PRINT iteration numbers 

fprintf (fiid, ' iterations = %d',itr); 

% CLOSE FILE "out.txt" 
f close (fiid) ; 

% FILE: obj ective_simult 

function L = obj ective_simult (a_inter) % FUNCTION DECLARE 

global pi ; 

global p2 ; 

global p3 ; 

global rmsl; 

global rms2 ; 

global rms3 ; 

global itr; 

ml=10; % Number of experimental tie lines 

m2=10; % Number of experimental tie lines 

m3=10; % Number of experimental tie lines 

% INCREASE Number of iterator by 1 

itr = itr+1 

% CONVERT the parameters in MATRIX form for System 1 
% a_interaction = [0 (1,2) (1,3) 

% (2,1) 0 (2,3) 

% (3,1) (3,2) 0 ] ; 

a_interactionl = [0 a_inter(l) a_inter(2); 

a_inter(3) 0 a_inter(4); 
a_inter(5) a_inter(6) 0] ; 

% CONVERT the parameters in MATRIX form for System 2 
a_interaction2 = [0 a_inter(7) a_inter(2); 

a_inter(8) 0 a_inter(9); 
a_inter(5) a_inter(10) 0] ,- 

% CONVERT the parameters in MATRIX form for System 3 
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a_interaction3 = [0 a_inter(ll) a_inter(2); 

a_inter(12) 0 a_inter (13 ) ; 
a_inter(5) a_inter(14) 0] ; 

% NOTE that parameters a_inter (2 ) and a_inter(5) are common 
% to all systems 

% OPEN FILE "exp_tiel.txt" for Experimental data 
fid = fopen ( ' exp_tiel . txt ' , ' r ' ) ; 

% READ tie line data in (6,m) MATRIX 
exptiel = f scanf (f id, ' %g %g',[6 inf] ) ; 

% CLOSE "exp_tiel.txt" 
f close (fid) ; 

% TAKE the transpose of MATRIX to make it (m,6) 
exp_tiel = exptiel ' ; 

% INITIATE sum's 

suml = 0 ; sum2 = 0 ; sum3 = 0 ; sum4 = 0 ; 

% INITIATE v's 
vl=0 ; v2=0 ; v3=0 ; v4=0 ; 

% CALL "estimation!" for tie line estimation 
% STORE them in a MATRIX 

est_tiel = estimation! (ml , exp_tiel , a_interactionl) ; 

% CALCULATE errors 

errorl = exp_tiel - est_tiel; 

% CALCULATE sum of errors 
for i = l:ml 

suml= suml+errorl (i , 1) ; 
sum2 = sum2 +error 1 ( i , 2 ) ; 
sum3= sumS+errorl (i , 4) ; 
sum4= sum4+errorl (i, 5) ; 

end 

% AVERAGE of errors 

suml = suml / ml ; sum2 = sum2 /ml ; sum3 = sum3 /ml ; 
sum4 = sum4/ml; 

% ESTIMATE variances 
for i = l:ml 

vl = vl+ (errorl (i , 1) -suml) ^2 ; 
v2 = v2+ (errorl (i, 2) -sum2) ^2 ; 
v3 = v3+ (errorl (i, 4) -sum3) ^2 ; 
v4 = v4+ (errorl (i , 5) -sum4) ^2 ; 

end 

vl=vl/ (ml-1) ;v2=v2/ (ml-1) ;v3=v3/ (ml-1) ;v4=v4/ (ml-1) ; 

% INITIATE 
tern = 0 ; 

% CALCULATE summation 
for i=l:ml 

tern = tern + (errorl (i , 1) ) ^2/vl+ (errorl (i , 2) ) ^2/v2+ .. . 

(errorl (i,4) ) '^2/v3+ (errorl (i, 5) ) ^2/v4 ; 
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end 


% CALCULATE the objective function 

LI = 2*ml*log (2*pi) + (ml/2) *log (vl*v2*v3*v4) + (1/2) *tem 
pi = LI; 
sum_square = 0 ; 

% CALCULATE the rmsd 
for i=l:ml 
for j=l:6 

sum_square = sum_square+ (est_tiel (i , j ) -exp_tiel (i, j ) ) ^2 ; 

end 

end 

rmsl = 100* (sum_square/ (6*ml) ) ^0 . 5 
%*************************************************** 

% Instructions are similar as above 
fid = fopen ( ' exp_tie2 . txt ' , ' r ' ) ; 
exptie2 = f scanf (f id, ' %g %g',[6 inf] ) ; 
f close (fid) ; 

exp_tie2 = exptie2 ' ; 

suml=0 ; sum2=0 ; sum3=0 ; sum4=0 ; 

vl = 0 ; v2 = 0 ; v3 = 0 ; v4 = 0 ; 

est_tie2 = estimation2 (m2 , exp_tie2 , a_interaction2 ) ; 
error2 = exp_tie2 - est_tie2; 
for i = l:m2 

suml= suml+error2 (i , 1) ; 
sum2= sum2+error2 (i,2) ; 
sum3= sum3+error2 (i , 4) ; 
sum4= sum4+error2 (i,5) ; 

end 

suml = suml/m2; sum2 = sum2/m2; sum3 = sum3/m2; 
sum4 = sum4 / m2 ; 
for i = l:m2 

vl = vl+ (error2 (i, 1) -suml) ^2 ; 
v2 = v2+ (error2 (i , 2 ) -sum2 ) ^2 ; 
v3 = v3+ (error2 (i, 4) -sum3) ^2 ; 
v4 = v4+ (error2 (i, 5) -sum4) ^2 ; 
end » 

vl=vl/ (m2-l) ;v2=v2/ (m2-l) ;v3=v3/ (m2-l) ;v4=v4/ (m2-l) ; 
tem = 0 ; 
for i=l:m2 

tem = tem + (error2 (i, 1) ) ^2/vl+ (error2 (i, 2) ) ^2/v2 . . . 

+ (error2 (i , 4 ) ) ^2/v3+ (error2 (i , 5) ) ^2/v4 ; 

end 

L2 = 2*m2*log (2*pi) + (m2/2) *log (vl*v2*v3*v4) + (1/2) *tem 
p2 = L2 ; 
sum_square = 0 ; 
for i=l:m2 
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for j=l:6 


sum_square = sum_square+ (est_tie2 (i, j ) - 
exp_tie2 ( i , j ) ) ''2 ; 
end 

end 

rms2 = 100* (sum_square/ (6*m2) ) ^0 . 5 

fid = fopen ( ' exp_tie3 . txt ' , ' r ' ) 

exptieS = f scanf (f id, ' %g %g',[6 inf ] ) ; 

f close (fid) ; 

exp_tie3 = exptie3 ' ; 

suml = 0 ; sum2 = 0 ; sum3 = 0 ; sum4 = 0 ; 

vl=0 ;v2=0 ;V3=0 ; v4=0 ; 

est_tie3 = estimation3 (m3 , exp_tie3 , a_interaction3 ) ; 
error3 = exp_tie3 - est_tie3; 
for i = l:m3 

suml= suml+error3 (i, 1) ; 
sum2 = sum2 +error3 ( i , 2 ) ; 
sum3 = sum3 +error 3 ( i , 4 ) ; 
sum4= sum4+error3 (i , 5) ; 

end 

suml = suml/m3; sum2 = sum2/m3; sum3 = sum3 /m3 ; 
sum4 = sum4/m3; 
for i = l:m3 

vl = vl+ (error! (i, 1) -suml) ^2 ; 
v2 = v2+ (error! (i, 2) -sum2) ^2 ; 
v3 = v3+ (error! (i , 4) -sum3) ^2 ; 
v4 = v4+ (error! (i , 5) -sum4) ^2 ; 

end 

vl=vl/ (m3-l) ;v2=v2/ (m3-l) ;v3=v3/ (m3-l) ;v4=v4/ (m3-l) ; 
tem = 0 ; 
for i=l:m3 

tem = tem + (error! (i , 1) ) ^2/vl+ (error! (i , 2) ) ^2/v2 .. . 

+ (error! (i,4) ) ^2/v!+ (error! (i, 5) ) ^2/v4; 

end 

L! = 2*m!*log (2*pi) + (m!/2) *log (vl*v2*v!*v4) + (1/2) *tem 
p! = L!; 
sum_square = 0 ; 
for i=l:m! 
for j=l:S 

sum_square = sum_square+ (est_tie! (i , j ) -exp_tie! (i , j ) ) ^2 ; 

end 

end 

rms! ~ 100* (sum__square/ (6*m! ) ) "^0 . 5 
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% SUM all the objective funations 
L = L1+L2+L3 


%************************************************* 

% FILE: estimation! .m 

^icicieicicicicicificicicicic-k-kic-k'kic-k'k'k-k'kicicicrkic'k'k-kic-k'k'kic-k-k'kicic-krkicik'kic 

% FUNCTION DECLARE 

function estimate = estimation! (m, exp_tie, a_interaction) 
% Rk for three components 
Rk = [4.4998 3.!905! 4.0356]; 

% Qk for three components 
Qk = [3.856 2.43.2]; 

% SPECIFY temperature 
temp =2 98.15 ; 
for i=!:m 

% ESTIMATE mid point of tie lines 

z_! = (exp_tie (i, !) +exp_tie (i,4) ) /2; 
z_2 = (exp_tie (i, 2) +exp_tie (i, 5) ) /2 ; 

% initial guess of a tie line 
XO = exp_tie (i , : ) ; 

% CALL "lle_equation" for setting LLE equations 
% CALL "fosolve" for Solving equations 
X = f solve (' lle_equation' , XO ,[],[], z_! , z_2 , Rk, Qk, 
a_int erection, temp) ; 

% STORE the tie lines in a MATRIX 
for j=!:6 

est (i, j ) = X(j ) ; 

end 

estimate = est; 

end 

% FILE: estimation2 . m 

% INSTRUCTION are similar to estimation! 

function estimate = estimation2 (m, exp_tie, a_interaction) 
Rk = [4.4998 3.92287 4.0356]; 

Qk = [3 . 856 2 . 968 3.2 ] ; 
temp =298.!5; 
for i=!:m 

z_! = (exp_tie (i, !) +exp_tie (i,4) ) /2 ; 
z_2 = (exp_tie (i, 2) +exp_tie (i, 5) ) /2 ; 

XO = exp_tie (i , : ) ; 
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X = f solve (' lle_equation' ,X0 ,[],[], z_l , z_2 , Rk, Qk, 

a_interaction, temp) ; 


double (X) ; 
for j=l:6 

est (i, j) = X(j) ; 

end 

estimate = est; 

end 

%**************-jlf**************'s*f*****************’^***** 

% FILE: estimations .m 

o 

% INSTRUCTION are similar to estimation! 

function estimate = estimations (m, exp_tie, a_interaction) 
Rk = [4.4998 4.65788 4.0S56]; 

Qk = [S.856 S.5S6 S.2 ] ; 
temp =298.15; 
for i=l:m 

z_l = (exp_tie (i , 1) +exp_tie (i , 4) ) /2 ; 
z_2 = (exp_tie (i , 2) +exp_tie (i , 5) ) /2 ; 

XO = exp_tie (i , : ) ; 

X = f solve (' lle_equation XO ,[],[], z_l , z_2 , Rk, Qk, 
a_int erection, temp) ; 

double (X) ; 
for j=l:6 

est (i , j ) = X ( j ) ; 

end 

estimate = est' ; 

end 
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A.4 Program to estimate temperature dependent parameters 




% FILE: initial_temp .m 
% INCLUDE objective_temp.m 
% INCLUDE estimation_temp .m 
% INCLUDE lle_equation.m 
% INCLUDE uniquac.m 
% INCLUDE unifac.m 
% INCLUDE nrtl.m 
% INCLUDE exptie.txt 


% same as above 
% same as above 
% same as above 
% same as above 
% It will contain 
temperature also 




global LK; 
global rms ; 
global itr; 

% START iteration 
itr = 0 ; 

% number of function evaluation 
options (14) = 5000; 

% INITIAL guess of parameters 
% a_inter0 = [a(l,2) a(l,3) b(l,2) b(l,3).. 

a(2,l) a(2,3) b(2,l) b(2,3) . . 
a(3,l) a(3,2) b(3,l) b(3,2) ] ; 
a_inter0 = [1083.86 -10.82902 246.3237 -2.25353 -1.1681 - 

11.3614 1.245044 -2.804212]; 

% CALL FILE: " act ivity_temp " for objective function 
% CALL "fmins" for minimization 

a_inter = fmins ( ' activity_temp ' , a_inter0 , options) 

% OPEN FILE "out.txt" for output 
fiid = fopen ( ' out . txt ' , ' a ' ) ; 

% PRINT initial guess 

fprintf (f iid, 'initial guess = %g %g %g %g %g %g %g %g %g %g %g 
%g\n' ,a_inter0) ; 

% PRINT parameters 

fprintf (fiid, 'parameters = %g %g %g %g %g %g %g %g %g %g %g 
%g\n' , a_inter) ; 

% PRINT objective functions 
fprintf (fiid, ' L = %g %g %g\n',LK); 

% PRINT rmsd's 

fprintf (fiid, 'pi = %g %g %g\n',rms); 

% PRINT no of iterations 

fprintf (fiid, ' iterations = %d',itr); 

% CLOSE FILE "out.txt" 
f close (fiid) ; 
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^ici(i(icici(icicieicicicicici(i(i<icici(icicici(icici(i(icic'kic'kic‘kic'kici(’k‘kicic‘k'k'k‘k-k'k-kicici(icicicicic 

% FILE: objective_teinp.m 

^******^*r***************r)«f**********'*^********‘********'*‘**‘*r**** 

function ob_fun = errorsq (a_inter) 
global LK; 
global rms; 
global itr; 

% SPECIFY number of tie lines for each system 
no_tie = [7,7,6]; 
no_tieline = no_tie; 

% SPECIFY number of temperatures for which data are 
no_data = 3 ; 
itr = itr+1 

% STORE the parameters in MATRIX form 
% a_interaction = [0 a{l,2)a(l,3)0 b(l,2)b(l,3); 

% a(2,l) 0 a(2,3)b(2,l) 0 b(2,3) ; 

% a(3,l)a(3,2) 0 b (3 , 1) b (3 , 2 ) 0 ]; 

a_interaction= [0 a_inter (1) a_inter (2) 0 a_inter (3 ) a_inter (4 ; 

a_inter{5)0 a_inter(6) a_inter(7) 0 a_inter(8); 
a_inter (9) a_inter (10) 0a_inter (11) a_inter (12) 0] 

% OPEN FILE "exptie.txt" 
fid = fopen ( ' exptie . txt ' , ' r ' ) ; 

% READ data in a MATRIX of (7,m) 
exptie = f scanf (f id, ' %g %g',[7 inf ] ) ,- 
% CLOSE FILE "out. txt" 
fclose (fid) ; 

% TRANSPOSE the MATRIX to make it (m,7) 
exp_tie = exptie'; 

% CALCULATE total number of tie lines 

term = 0 ; 

for k=l : no_data 

term = term + no_tie (k) ; 
no_tie (k) = term; 

end 

% CALL "estimation_temp" for tie line estimation 
est_tie=estimation (a_interaction, exp_tie, no_tie (no_data) ) ; 
% STORE experimental tie line data in MATRIX 
exp_tie = exp_tie ( : , 1 : 6) ; 

% ESTIMATE the errors 

error = exp_tie - est_tie; 

no_tie ; 

terml = 0 ; 

term2 = 0 ; 

term3 = 0 ; 

objective = 0; 
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% ESTIMATE all objective functions for all temperature data 
for k = l:no_data 

sum_square (k) = 0; 
suml = 0 ; sum2 = 0 ; sum3 = 0 ; sum4 = 0 ; 
vl=0 ; v2=0 ; v3=0 ; v4=0 ; 
m = no_tieline (k) ; 
for i = terml+1 :no_tie (k) 
terml = no_tie (k) ; 
suml= suml+error (i , 1) ; 
sum2 = sum2 +error ( i , 2 ) ; 
sum3 = sum3 +error ( i , 4 ) ; 
sum4= sum4+error (i , 5) ; 

end 

suml = suml/m; sum2 = sum2/m; 
sum3 = sum3/m; sum4 = sum4/m; 
for i = term2+l :no_tie (k) 
term2=no_tie (k) ; 
vl = vl+ (error (i , 1) -suml) ^2 ; 
v2 = v2+ (error (i,2) -sum2) ^2; 
v3 = v3+ (error (i , 4) -sum3 ) "^2 ; 
v4 = v4+ (error (i , 5) -sum4) ^2 ; 

end 

vl=vl/ (m-1) ;v2=v2/ (m-1) ;v3=v3/ (m-1) ;v4=v4/(m-l) ; 
tem = 0 ; 

for i=term3+l : no_tie (k) 
term3=no_tie (k) ; 

tem = tem +(error(i,l))'^2/vl+(error(i,2))^2/v2 + 

(error (i , 4) ) ^2/v3 + (error (i , 5) ) ^2/v4 ; 

for j =1:6 

sum_square (k) = sum_square (k) + (error (i , j ) ) ^2 ; 

end 


end 

L(k) = 2*m*log (2*pi) + (m/2) *log (vl*v2*v3*v4) + (1/2) *tem; 

objective = objective + L(k); 

rmsd(k) = 100* (sum_square (k) / (m*6) ) ^0 . 5 ; 

LK(k) = L (k) ; 
rms (k) =rmsd (k) ; 

end 

L 

% ob_fun is sum of all objective functions 

ob_fun = objective 

rmsd 
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% FILE: estimaton_temp . m 

%****************************************************** 
function estimate = estimation (a_interaction, exp_tie,m) 
for i=l:m 

z_l = (exp_tie (i , 1) +exp_tie (i, 4) ) /2 ; 
z_2 = (exp_tie (i, 2) +exp_tie (i, 5) ) /2 ; 

XO = exp_tie (i , 1 : 6) ; 
temp = exp_tie ( i , 7 ) ; 

X = fsolve (' fun' ,X0, [],[], z_l, z_2,a_interaction, temp) ; 
double (X) ; 
for j=l:6 

est (i, j ) = X(j ) ; 

end 

estimate = est; 
end 
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